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Abstract. The connection between the standard i?°°-problcm in control the- 
ory and Nevanlinna-Pick interpolation in operator theory was established in 
the 1980s, and has led to a fruitful cross-pollination between the two fields 
since. In the meantime, research in i?°°-control theory has moved on to the 
study of robust control for systems with structured uncertainties and to var- 
ious types of multidimensional systems, while Nevanlinna-Pick interpolation 
theory has moved on independently to a variety of multivariable settings. Here 
we review these developments and indicate the precise connections which sur- 
vive in the more general multidimensional/multivariable incarnations of the 
two theories. 



1. Introduction 

Starting in the early 1980s with the seminal paper [139] of George Zanies, there 
occurred an active interaction between operator theorists and control engineers in 
the development of the early stages of the emerging theory of iJ°°-control. The 
cornerstone for this interaction was the early recognition by Francis-Hclton-Zames 
[65] that the simplest case of the central problem of i7°°-control (the sensitivity 
minimization problem) is one and the same as a Nevanlinna-Pick interpolation 
problem which had already been solved in the early part of the twentieth century 
(see |1101 llOSj ). For the standard problem of i7°°-control it was known early on 
that it could be brought to the so-called Model-Matching form (see [131 IM])- In 
the simplest cases, the Model-Matching problem converts easily to a Nevanlinna- 
Pick interpolation problem of classical type. Handling the more general problems 
of _ff °°-control required extensions of the theory of Nevanlinna-Pick interpolation 
to tangential (or directional) interpolation conditions for matrix-valued functions; 
such extensions of the interpolation theory were pursued by both engineers and 
mathematicians (see e.g. [23 liHl ISOl IHEl [HI] ) ■ Alternatively, the Model-Matching 
problem can be viewed as a Sarason problem which is suitable for application of 
Commutant Lifting theory (see [125, ,62 ). The approach of ^64j used an additional 
conversion to a Nehari problem where existing results on the solution of the Nehari 
problem in state-space coordinates were applicable (see |69l [33]). The book of 
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Francis [SJ was the first book on _ff ""-control and provides a good summary of tlie 
state of the subject in 1987. 

While there was a lot of work emphasizing the connection of the iJ°°-problem 
with interpolation and the related approach through J-spectral factorization f [26l 
[901 [911 [Ml ill [331 [24]), we should point out that the final form of the 
theory parted ways with the connection with Nevanlinna-Pick interpolation. When 
calculations were carried out in state-space coordinates, the reduction to Model- 
Matching form via the Youla-Kucera parametrization of stabilizing controllers led 
to inflation of state-space dimension; elimination of non-minimal state-space nodes 
by finding pole-zero cancellations demanded tedious brute-force calculations (see 
[nniini])- a direct solution in state-space coordinates (without reduction to Model- 
Matching form and any explicit connection with Nevanlinna-Pick interpolation) was 
finally obtained by Ball-Cohen ^24J (via a J-spectral factorization approach) and 
in the more definitive coupled-Riccati-equation form of Doyle-Glover-Khargonekar- 
Francis |54) . This latter paper emphasizes the parallels with older control paradigms 
(e.g., the Linear-Quadratic-Gaussian and Linear-Quadratic- Regulator problems) 
and obtained parallel formulas for the related iJ-^-problem. The J-spectral fac- 
torization approach was further developed in the work of Kimura, Green, Glover, 
Limebeer, and Doyle [ST llTDllTTj . A good review of the state of the theory to this 
point can be found in the books of Zhou-Doyle-Glover |141j and Green-Limebeer 

The coupled-Riccati-equation solution however has now been superseded by the 
Linear-Matrix-Inequality (LMI) solution which came shortly thereafter; we mention 
specifically the papers of Iwasaki-Skelton [7S] and Gahinet-Apkarian This so- 
lution does not require any boundary rank conditions entailed in all the earlier 
approaches and generalizes in a straightforward way to more general settings (to be 
discussed in more detail below) . The LMI form of the solution is particularly appeal- 
ing from a computational point of view due to the recent advances in semidefinite 
programming (see [68 ). The book of DuUerud-Paganini [57] gives an up-to-date 
account of these latest developments. 

Research in iJ°°-control has moved on in a number of different new directions, 
e.g., extensions of the iJ°°-paradigm to sampled-data systems [33, nonlinear sys- 
tems |126j . hybrid systems [23], stochastic systems |76], quantum stochastic sys- 
tems [79], linear repetitive processes [123 , as well as behavioral frameworks jl34| . 
Our focus here will be on the extensions to robust control for systems with struc- 
tured uncertainties and related i/°°-control problems for multidimensional (A'^- 
D) systems — both frequency-domain and state-space settings. In the meantime, 
Nevanlinna-Pick interpolation theory has moved on to a variety of multivariable 
settings (polydisk, ball, noncommutative polydisk/ball); we mention in particular 
the papers [Il[49l[ri3, 3, 35, 19, 20, 21, 22,^. 

As the transfer function for a multidimensional system is a function of several 
variables, one would expect that the same connections familiar from the 1-D/single- 
variable case should also occur in these more general settings; however, while there 
had been some interaction between control theory and several-variable complex 
function theory in the older area of systems over rings (see [83l [85l to this 

point, with a few exceptions |731 1741 132j . there has not been such an interaction 
in connection with iJ°°-control for A^-D systems and related such topics. With 
this paper we wish to make precise the interconnections which do exist between 
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the _ff°°-theory and the interpolation theory in these more general settings. As we 
shall see, some aspects which are taken for granted in the 1-D/single- variable case 
become much more subtle in the A^-D/multivariable case. Along the way we shall 
encounter a variety of topics that have gained attention recently, and sometimes 
less recently, in the engineering literature. 

Besides the present Introduction, the paper consists of five sections which we 
now describe: 

(1) In Section [5] we lay out four specific results for the classical 1-D case; these 
serve as models for the type of results which we wish to generalize to the A^- 
D/multivariable settings. 

(2) In Section [3] we survey the recent results of Quadrat [TTZl EIHl El \IM 
11211 1122] on internal stabilization and parametrization of stabilizing controllers in 
an abstract ring setting. The main point here is that it is possible to parametrize 
the set of all stabilizing controllers in terms of a given stabilizing controller even 
in settings where the given plant may not have a double coprime factorization — 
resolving some issues left open in the book of Vidyasagar ^136j . In the case where 
a double-coprime factorization is available, the parametrization formula is more 
efficient. Our modest new contribution here is to extend the ideas to the setting 
of the standard problem of i/°°-control (in the sense of the book of Francis [B11 ) 
where the given plant is assumed to have distinct disturbance and control inputs 
and distinct error and measurement outputs. 

(3) In Section m we look at the internal-stabilization/i7°°-control problem for 
multidimensional systems. These problems have been studied in a purely frequency- 
domain framework (see [HllIM]) as well as in a state-space framework (see [STl [551 
[5B ). In Subsection 14. 11 we give the frequency-domain formulation of the problem. 
When one takes the stable plants to consist of the ring of structurally stable ra- 
tional matrix functions, the general results of Quadrat apply. In particular, for 
this setting stabilizability of a given plant implies the existence of a double co- 
prime factorization (see |119j ). Application of the Youla-Kucera parametrization 
then leads to a Model-Matching form and, in the presence of some boundary rank 
conditions, the i7°°-problem converts to a polydisk version of the Nevanlinna-Pick 
interpolation problem. Unlike the situation in the classical single-variable case, 
this interpolation problem has no practical necessary-and-sufficient solution crite- 
rion and in practice one is satisfied with necessary and sufficient conditions for the 
existence of a solution in the more restrictive Schur-Agler class (see [Tl \S[ [55]). 

In Subsection l4.2l we formulate the internal-stabilization/77°°-control problem in 
Givone-Roesser state-space coordinates. We indicate the various subtleties involved 
in implementing the state-space version [104|, [85] of the double-coprime factorization 
and associated Youla-Kucera parametrization of the set of stabilizing controllers. 
With regard to the i/°°-control problem, unlike the situation in the classical 1-D 
case, there is no useable necessary and sufficient analysis for solution of the problem; 
instead what is done (see e.g. [55l[56]) is the use of an LMI/Bounded- Real-Lemma 
analysis which provides a convenient set of sufficient conditions for solution of the 
problem. This sufficiency analysis in turn amounts to an N-D extension of the LMI 
solution ^78| 166] of the 1-D i/°°-control problem and can be viewed as a necessary 
and sufficient analysis of a compromise problem (the "scaled" i7°°-problem). 

While stabilization and i?°°-control problems have been studied in the state- 
space setting HlJ [551 [56] and in the frequency-domain setting [92l [93] separately. 
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there does not seem to have been much work on the precise connections between 
these two settings. The main point of Subsection 14.31 is to study this relationship; 
while solving the state-space problem implies a solution of the frequency-domain 
problem, the reverse direction is more subtle and it seems that only partial re- 
sults are known. Here we introduce a notion of modal stabilizability and modal 
detectability (a modification of the notions of modal controllability and modal ob- 
servability introduced by Kung-Levy-Morf-Kailath ^88j ) to obtain a partial result on 
relating a solution of the frequency-domain problem to a solution of the associated 
state-space problem. This result suffers from the same weakness as a correspond- 
ing result in \88\ : just as the authors in [88j were unable to prove that minimal 
(i.e., simultaneously modally controllable and modally observable) realizations for 
a given transfer matrix exist, so also we are unable to prove that a simultaneously 
modally stabilizable and modally detectable realization exists. A basic difficulty in 
translating from frequency-domain to state-space coordinates is the failure of the 
State-Space-Similarity theorem and related Kalman state-space reduction for A'^-D 
systems. Nevertheless, the result is a natural analogue of the corresponding 1-D 
result. 

There is a parallel between the control-theory side and the interpolation-theory 
side in that in both cases one is forced to be satisfied with a compromise solution: 
the scaled- problem on the control-theory side, and the Schur-Agler class (rather 
than the Schur class) on the interpolation-theory side. We include some discussion 
on the extent to which these compromises are equivalent. 

(4) In Section [5] we discuss several 1-D variations on the internal-stabilization 
and /J°°-control problem which lead to versions of the A^-D/multivariable prob- 
lems discussed in Section S) It was observed early on that an i7°°-controller has 
good robustness properties, i.e., an i7°°-controller not only provides stability of the 
closed-loop system associated with the given (or nominal) plant for which the con- 
trol was designed, but also for a whole neighborhood of plants around the nominal 
plant. This idea was refined in a number of directions, e.g., robustness with respect 
to additive or multiplicative plant uncertainty, or with respect to uncertainty in 
a normalized coprime factorization of the plant (see |100j ). Another model for an 
uncertainty structure is the Linear-Fractional- Transformation (LFT) model used 
by Doyle and coworkers (see [97l [98]). Here a key concept is the notion of struc- 
tured singular value fJ-{A) for a finite square matrix A introduced by Doyle and 
Safonov [521 1124] which simultaneously generalizes the norm and the spectral ra- 
dius depending on the choice of uncertainty structure (a C*-algebra of matrices 
with a prescribed block-diagonal structure); we refer to |107j for a comprehensive 
survey. If one assumes that the controller has on-line access to the uncertainty 
parameters one is led to a gain-scheduling problem which can be identified as the 
type of multidimensional control problem discussed in Section 14.21 — see [1061 [T8] ; 
we survey this material in Subsection 15.11 In Subsection 15.21 we review the purely 
frequency-domain approach of Helton |73l[74] toward gain-scheduling which leads to 
the frequency-domain internal-stabilization/if °°-control problem discussed in Sec- 
tion [4Tj Finally, in Section 15.31 we discuss a hybrid frequency-domain/state-space 
model for structured uncertainty which leads to a generalization of Nevanlinna-Pick 
interpolation for single-variable functions where the constraint that the norm be 
uniformly bounded by 1 is replaced by the constraint that the /z-singular value be 
uniformly bounded by 1 ; this approach has only been analyzed for very special cases 
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of the control problem but does lead to interesting new results for operator theory 
and complex geometry in the work of Bercovici-Foias-Tannenbaum [3H1 [Ml IIHI SI] , 
Agler- Young [SlilHIHlilliniinillllll^, Huang-Marcantognini- Young f77], and 
Popescu |114j . 

(5) The final Section [6] discusses an enhancement of the LFT-model for struc- 
tured uncertainty to allow dynamic time-varying uncertainties. If the controller is 
allowed to have on-line access to these more general uncertainties, then the solution 
of the internal-stabilization/i/'^-control problem has a form completely analogous 
to the classical 1-D case. Roughly, this result corresponds to the fact that, with this 
noncommutative enhanced uncertainty structure, the a priori upper bound /i(A) 
for the structured singular value /i(A) is actually equal to /i(A), despite the fact 
that for non-enhanced structures, the gap between /i and Jl can be arbitrarily large 
(see |133p . In this precise form, the result appears for the first time in the the- 
sis of Paganini |108j but various versions of this type of result have also appeared 
elsewhere (see [37l |42l [60l [99l I129| ). We discuss this enhanced noncommutative 
LFT-model in Subsection 16.11 In Subsection 16.21 we introduce a noncommutative 
frequency-domain control problem in the spirit of Chapter 4 of the thesis of Lu 
[96j . where the underlying polydisk occurring in Section [4. II is now replaced by the 
noncommutative polydisk consisting of all d-tuples of contraction operators on a 
fixed separable infinite-dimensional Hilbert space K. and the space of i/°°-functions 
is replaced by the space of scalar multiples of the noncommutative Schur- Agler class 
introduced in |28j . Via an adaptation of the Youla-Kucera parametrization of stabi- 
lizing controllers, the internal-stabilization/77°°-control problem can be reduced to 
a Model-Matching form which has the interpretation as a noncommutative Sarason 
interpolation problem. In the final Subsection 16.31 we show how the noncommu- 
tative state-space problem is exactly equivalent to the noncommutative frequency- 
domain problem and thereby obtain an analogue of the classical case which is much 
more complete than for the commutative-variable case given in Section 14.31 In 
particular, if the problem data are given in terms of state-space coordinates, the 
noncommutative Sarason problem can be solved as an application of the LMI so- 
lution of the i7°°-problem. While there has been quite a bit of recent activity on 
this kind of noncommutative function theory (see e.g. [l4l[22ll75lf82l[TT5llll6p . the 
noncommutative Sarason problem has to this point escaped attention; in particu- 
lar, it is not clear how the noncommutative Nevanlinna-Pick interpolation problem 
studied in [22j is connected with the noncommutative Sarason problem. 

Finally we mention that each section ends with a "Notes" subsection which dis- 
cusses more specialized points and makes some additional connections with existing 
literature. 

Acknowledgement. The authors thank Quanlei Fang and Gilbert Groenewald for 
the useful discussions in an early stage of preparation of the present paper. We also 
thank the two anonymous reviewers for their thorough readings of the first version 
and constructive suggestions for the preparation of the final version of this paper. 

2. The 1-D systems/single- variable case 

Let C[z] be the space of polynomials with complex coefficients and C(z) the 
quotient field consisting of rational functions in the variable z. Let TZH°° be the 
subring of stable elements of C(z) consisting of those rational functions which are 
analytic and bounded on the unit disk D, i.e., with no poles in the closed unit disk 
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D. We assume to be given a plant G = [021022]-^®^"*'^®^ which is 
given as a block matrix of appropriate size with entries from C(z). Here the spaces 
W, Z and y have the interpretation of control- signal space, disturbance- signal 
space, error-signal space and measurement- signal space, respectively, and consist of 
column vectors of given sizes nu, jt-vVj and ny, respectively, with entries from 
C(z). For this plant G we seek to design a controller K: y —> U, also given as a 
matrix over C(z), that stabilizes the feedback system Y,{G,K) obtained from the 
signal- flow diagram in Figure [T] in a sense to be defined precisely below. 
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Figure 1. Feedback with tap signals 



Note that the various matrix entries Gy of G are themselves matrices with entries 
from C(z) of compatible sizes (e.g., Gn has size nz x riw) and X is a matrix over 
C{z) of size ni( x ny. 

The system equations associated with the signal-flow diagram of Figure [T] can 
be written as 

7 -G12 

/ -~K u ^ I vi . (2.1) 
-G22 / 

Here vi and V2 are tap signals used to detect stability properties of the internal 
signals u and y. We say that the system S(G, if) is well-posed if there is a well- 

defined map from "i to « . It follows from a standard Schur complement 
L "2 J L f J 

computation that the system is well-posed if and only if dct(/ — G22K) 7^ 0, and 
that in that case the map from 
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We say that the system J^{G,K) is internally stable if Y.{G,K) is well-posed and, 
in addition, if the map e(G, K) maps TZH^ © UH^ © TZH^ into TZH^ © TZH^ © 
TZHy' , i.e., stable inputs w,fi,W2 are mapped to stable outputs z,u,y. Note that 
this is the same as the condition that the entries of S(G, K) be in TZH°° . 

We say that the system S(G', K) has performance if S(G, is internally stable 
and in addition the transfer function Tzw from w to z has supremum-norm over the 
unit disk bounded by some tolerance which we normalize to be equal to 1: 

II '^xii? 1 1 oo • — 

sup{||r,^(A)||: AgB}<1. 

Here ||T'2^(A)|| refers to the induced operator norm, i.e., the largest singular value 
for the matrix T^wW- We say that the system S(G, K) has strict performance if in 
addition HTz^Hoo < 1- The stabilization problem then is to describe all (if any exist) 
internally stabilizing controllers K for the given plant G, i.e., all K £ C(z)""^"^ so 
that the associated closed-loop system I]{G,K) is internally stable. The standard 
H°° -problem is to find all internally stabilizing controllers which in addition achieve 
performance HT^^jIIcx) < 1- The strictly suboptinial i/°°-problem is to describe all 
internally stabilizing controllers which also achieve strict performance HTz^, ||oo < 1. 

2.1. The model-matching problem. Let us now consider the special case where 
G22 — 0, so that G has the form G = [ '"^J^ ] . In this case well-posedness is 
automatic and 8(G, /f) simplifies to 



e(G,X) = 



Gil + G12KG21 G12 G12K 
KG21 I K 

G21 / 



Thus internal stability for the closed- loop system S(G, K) is equivalent to stability 
of the four transfer matrices Gn, G12, G21 and K. Hence internal stabilizability of 
G is equivalent to stability of Gn, G12 and G21; when the latter holds a given K 
internally stabilizes G if and only if K itself is stable. 

Now assume that Gn, G12 and G21 are stable. Then the i/°°-performance prob- 
lem for G consists of finding stable K so that ||Gii -I- Gi2iirG2i ||oo < 1- Following 
the terminology of [64], the problem is called the Model-Matching Problem. Due to 
the influence of the paper [125j . this problem is usually referred to as the Sarason 
problem in the operator theory community; in ^125j it is shown explicitly how the 
problem can be reduced to an interpolation problem. 

In general control problems the assumption that G22 = is an unnatural as- 
sumption. However, after making a change of coordinates using the Youla-Kucera 
parametrization or the Quadrat parametrization, discussed below, it turns out that 
the general _ff°°-problem can be reduced to a model-matching problem. 

2.2. The frequency-domain stabilization and H°° problem. The following 
result on characterization of stabilizing controllers is well known (sec e.g. [64 or 
[136[ I137j for a more general setting) . 

Theorem 2.1. Suppose that we are given a rational matrix function G = [q^^ ] 
of size {uz + ny) x (nw + nu) with entries in C(z) as above. Assume that G is 
stabilizable, i.e., there exists a rational matrix function K of size nu x ny so that 
the nine transfer functions in p.2p are all stable. Then a given rational matrix 
function K stabilizes G if and only if K stabilizes G22, i-e., Q{G,K) in (|2.2p is 
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stable if and only if 

err r.^ . + k{i -G22K)-^ 
{1-KG22)-' 

G22{I - KG22)-^ I + G22{I - KG22) 

is stable. Moreover, if we are given a double coprime factorization for G22, i-e., 
stable transfer matrices D, N , X , Y, D, N , X and Y so that the determinants of 
D, D, X and X are all nonzero (in TZH°° ) and 
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(2.3) 



(such double coprime factorizations always exists since TZH°° is a Principal Ideal 
Domain), then the set of all stabilizing controllers K is given by either of the for- 
mulas 

K={Y + DA){X + NA)-^ = iX + AN^^iY + AD), 
where A is a free stable parameter from T^H^fjj y-^ such that det(X + NA) ^ or 
equivalently det{X + AN) ^ 0. 

Through the characterization of the stabihzing controUers, those controllers that, 
in addition, achieve performance can be obtained from the solutions of a Model- 
Matching/Sarason interpolation problem. 

Theorem 2.2. Assume that G £ C(z)("^+"^)''(""'+"") is stabilizable and that G22 
admits a double coprime factorization (|3.9p . Let K € C(z)""^"y. Then K is a 
solution to the standard H°° problem for G if and only if 



{Y + DA){X + 7VA)~i = (X + AN)^^{Y + AD), 
^c{u y) ^^'^^ det(X + NA) ^ 0, or equivalently det(X + AiV) ^ 0, 



K 

where A G TZH^ 

is any solution to the Model-Matching/Sarason interpolation problem for Gu, G12 
and G21 defined by 

Gn :— Gil 

i.e., so that 



G12YDG21, Gi2:—Gi2D, G21 DG21, 
IIG11 + G12AG21IU < 1. 



We note that in case G12 is injective and G21 is surjective on the unit cir- 
cle, by absorbing outer factors into the free parameter A we may assume without 
loss of generality that G12 is inner (i.e., Gi2{z) is isometric for z on unit cir- 
cle) and G21 is co-inner (i.e., G2i{z) is coisometric for z on the unit circle). Let 
T: Ly^ Q G2iHy^ L\ G^H^ be the compression of multiplication by Gu to 
the spaces G^ii?^^ and L% Gi2H^, i.e., T = Pl^^qg^^h^.^uIli^qg-hI^- 
Then, as a consequence of the Commutant Lifting theorem (see [531 Corollary 10.2 
pages 40-41]), one can see that the strict Model-Matching/Sarason interpolation 
problem posed in Theorem 12.21 has a solution if and only if ||r|lop < 1. Alter- 
natively, in case G12 and G21 are square and invertible on the unit circle, one 
can convert this Model-Matching/Commutant-Lifting problem to a bitangential 



Nevanlinna-Pick interpolation problem (see Theorem 16.9.3]), a direct general- 
ization of the connection between a model-matching/Sarason interpolation problem 
with Nevanlinna-Pick interpolation as given in [1251 165 ] for the scalar case, but we 
will not go into the details of this here. 

2.3. The state-space approach. We now restrict the classes of admissible plants 
and controllers to the transfer matrices whose entries are in C(2;)o, the space of 
rational functions without a pole at (i.e., analytic in a neighborhood of 0). In 
that transfer matrix F :IA ^ y with entries in C(z)o admits a state-space 

realization: There exists a quadruple {A, B, C, D} consisting of matrices whose sizes 
are given by 
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■ X ' 
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(2.4) 



where the state-space X is finite dimensional, so that 



for z in a neighborhood of 0. Sometimes we consider quadruples {A, B, C, D} of 
operators, of compatible size as above, without any explicit connection to a transfer 
matrix, in which case we just speak of a realization. 

Associated with the realization {A, B, C, D} is the linear discrete-time system 
of equations 



The system S and function F are related through the fact that F is the transfer- 
function of S. The two-by-two matrix (|2.4[) is called the system matrix of the system 



For the rest of this section we shall say that an operator A on a finite-dimensional 
state space X is stable if all its eigenvalues are in the open unit disk, or, equivalcntly, 
]|yl"a;]| — > as n — > oo for each x € X. The following result deals with two key 
notions for the stabilizability problem on the state-space level. 

Theorem 2.3. (I) Suppose that {A, B} is an input pair, i.e., A, B are operators 
with A: X X and B: lA ^ X for a finite- dimensional state space X and a 
finite- dimensional input space U. Then the following are equivalent: 

(1) {A, B} is operator-stabilizable, i.e., there exists a state-feedback operator 
F : X —> lA so that the operator A + BF is stable. 

(2) {A,B} «s Hautus-stabilizable, i.e., the matrix pencil [l — zA _B] is sur- 
jective for each z in the closed unit disk D. 

(3) The Stein inequality 



has a positive- definite solution X . Here F < for a square matrix F means 
that — F is positive definite. 

(II) Dually, if {C,A} is an output pair, i.e., C,A are operators with A: X ^ X 
and C : X y for a finite- dimensional state space X and a finite- dimensional 
output space y, then the following are equivalent: 

(1) {C, ^} is operator-detectable, i.e., there exists an output-injection operator 
L : y ^ X so that A -\- LC is stable. 



F{z) ^ D + zC{I - zAy^B 




Ax{n) + Bu{n), 
Cx{n) + Du{n). 



[n e Z+) 



AX A* - X - BB* < 
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is injective for 



(2) {C, ^} is Hautus-detectable, i.e., the matrix pencil J"*] 
all z in the closed disk D. 

(3) The Stein inequality 

A*YA-Y - C*C <0 
has a positive definite solution Y . 

When the input pair {A, B} satisfies any one (and hence all) of the three equiv- 
alent conditions in part (I) of Theorem 12. 3[ we shall say simply that {A, B} is 
stabilizable. Similarly, if (C, A) satisfies any one of the three equivalent conditions 
in part (II), we shall say simply that {C, A} is detectable. Given a realization 
{A, B, C, D}, we shall say that {A, B, C, D} is stabilizable and detectable if {A, B} 
is stabilizable and {C, A} is detectable. 

In the state-space formulation of the internal stabilization/// ""-control problem, 
one assumes to be given a state-space realization for the plant G: 

D21 D22 



G{z) 



+ z 



C2 



{I^zA)-^[B^ B2] 



(2.5) 



where the system matrix has the form 

A Bi B2 



Ci 

C2 



D21 



D12 
D22 





' X' 




'X 




w 




z 




u 




y 



(2.6) 



One then seeks a controller K which is also given in terms of a state-space realization 
K{z) = Dk + zCk{I - zAk)-^Bk 

which provides internal stability (in the state-space sense to de defined below) 
and/or iJ°°-pcrformance for the closed-loop system. Well-posedness of the closed- 
loop system is equivalent to invertibility of / — D22Dk. To keep various formulas 
affine in the design parameters Ak, Bk, Ck, Dk, it is natural to assume that D22 = 
0; this is considered not unduly restrictive since under the assumption of well- 
posedness this can always be arranged via a change of variables (see [78]). Then 
the closed loop system &{G, K) admits a state space realization {A^, B^i, Cd, Dd} 
given by its system matrix 



Aci 
Gel 



Bel 
Del 



A- 


f B2DKC2 


B2CK 


Bi + B2DKD21 




BkG2 


Ak 


BkD21 


. ^1- 


f D12DKG2 


D12GK 


D11+D12DKD21 



(2.7) 



B2DKC2 

BkC2 



and internal stability (in the state-space sense) is taken to mean that Ad = 
should be stable, i.e., all eigenvalues are in the open unit disk. 

The following result characterizes when a given G is internally stabilizable in the 
state-space sense. 

Theorem 2.4. (See Proposition 5.2 in |57j . ) Suppose that we are given a system 
matrix as in (|2.6p with D22 = with associated transfer matrix G as in (j2.5p . Then 
there exists a K{z) = Dk+zGk{I^zAk)^^Bk which internally stabilizes G (in the 
state-spaces sense) if and only if {A, B2} is stabilizable and {G2,A} is detectable. 
In this case one such controller is given by the realization {Afc, Bk,Gk, T>k} with 
system matrix 



B2 Cn 
Ak 



' Ak 


Bk ' 




'A + B2F + LG2 -L 




Dk . 




F 
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where F and L are state-feedback and output-injection operators chosen so that 
A + B2F and A + LC2 are stable. 



In addition to the state-space version of the stabihzabihty problem we also con- 
sider a (strict) state-space H°° problem, namely to find a controller K given by a 
state-space realization {Ak, Bk,Ck, Dk} of compatible size so that the transfer- 
function Tzui of the closed loop system, given by the system matrix (|2.7|) , is stable 
(in the state-space sense) and has a supremum norm ||r2u,||oo of at most 1 (less 
than 1). 

The definitive solution of the _ff °°-control problem in state-space coordinates for 
a time was the coupled-Riccati-equation solution due to Doyle-Glover-Khargonekar- 
Francis [53] . This solution has now been superseded by the LMI solution of Gahinet- 
Apkarian [55] which can be stated as follows. Note that the problem can be solved 
directly without first processing the data to the Model-Matching form. 

Theorem 2.5. Let {A, B,C, D} ^ {A,[B^ b^] ,[c^ C2] , [ ]} be a given re- 

alization. Then there exists a solution for the strict state-space H°° -control problem 
associated with {A, B, C, D} if and only if there exist positive-definite matrices X, Y 
satisfying the LMIs 





No 




AY A* - Y 
C'lYA* 
Bl 

'A*XA - X 
BIX A 



AYCl 

CiYCl - 

n* 
^11 

A*XBi 
BlXBi - 



I 



-I 
D* 



and the coupling condition 



X I 
I Y 



> 0. 





No 




< 0, y > 0, 



< 0, x>o, 



(2.8) 
(2.9) 

(2.10) 



Here Nc and No are matrices chosen so that 

Nc is injective and Im Nc — Ker [-B2 £'12] '^'^'^ 
No is injective and ImiVo = Ker [C2 D2i\ . 



We shall discuss the proof of Theorem 1 2. 5 1 in Section [3?2l below in the context of 
a more general multidimensional-system 77°°-control problem. 

The next result is the key to transferring from the frequency-domain version of 
the internal-stabilization/i7°°-control problem to the state-space version. 

Theorem 2.6. (See Lemma 5.5 in |57| .) Suppose that the realization {A, B2, C2, 0} 

for the plant G22 o,nd the realization {Ak, Bk,Ck, L)k} for the controller K are 
both stabilizable and detectable. Then K internally stabilizes G22 in the state-space 
sense if and only if K stabilizes G22 in the frequency- domain sense, i.e., the closed- 
loop matrix Ad = [ ^^Bk<!2'^^ ^A^"^ ] ^•^ stable if and only if the associated transfer 
matrix 



e{G22,K) = 



I Dk 
/ 



DkC2 Gk 
G2 



(/ - zAoiY 



B2 B2DK 
Bk 



has all matrix entries in TZH^ 
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2.4. Notes. In the context of the discussion immediately after the statement of 
Theorem 12.21 in case G12 and/or G21 drop rank at points on the unit circle, the 
Model-Matching problem in Theorem 12.21 may convert to a boundary Nevanlinna- 
Pick interpolation problem for which there is an elaborate specialized theory (see 
e.g. Chapter 21 of [26 and the more recent [43]). However, if one sticks with the 
strictly suboptimal version of the problem, one can solve the problem with the 
boundary interpolation conditions if and only if one can solve the problem without 
the boundary interpolation conditions, i.e., boundary interpolation conditions are 
irrelevant as far as existence criteria are concerned. This is the route taken in the 
LMI solution of the iJ°°-problem and provides one explanation for the disappear- 
ance of any rank conditions in the formulation of the solution of the problem. For 
a complete analysis of the relation between the coupled-Riccati-equation of [SI] 
versus the LMI solution of [55], we refer to |127j . 

3. The fractional representation approach to stabiliz ability and 

performance 

In this section we work in the general framework of the fractional representa- 
tion approach to stabilization of linear systems as introduced originally by Desoer, 
Vidyasagar and coauthors [SO] I137j in the 1980s and refined only recently in the 
work of Quadrat [1181 11211 1122] . For an overview of the more recent developments 
we recommend the survey article |117) and for a completely elementary account of 
the generalized Youla-Kucera parametrization with all the algebro-geometric inter- 
pretations stripped out we recommend |120j . 

The set of stable single-input single-output (SISO) transfer functions is assumed 
to be given by a general ring A in place of the ring TZH°° used for the classical 
case as discussed in Section (2) the only assumption which we shall impose on A 
is that it be a commutative integral domain. It therefore has a quotient field 
K := (5(A) = {n/d: d,n G A,d ^ 0} which shall be considered as the set of all 
possible SISO transfer functions (or plants). Examples of A which come up include 
the ring Rs(2) of real rational functions of the complex variable z with no poles in 
the closed right half plane, the Banach algebra RH°°{C+) of all bounded analytic 
functions on the right half plane C4- which are real on the positive real axis, and 
their discrete-time analogues: (1) real rational functions with no poles in the closed 
unit disk (or closed exterior of the unit disk depending on how one sets conventions) , 
and (2) the Banach algebra RH°°{D) of all bounded holomorphic functions on the 
unit disk D with real values on the real interval (—1, 1). There are also Banach 
subalgebras of RH°°{C+) or RH°°{ISi) (e.g., the Wiener algebra and its relatives 
such as the Callier-Desoer class — see [48]) which are of interest. In addition to these 
examples there are multivariable analogues, some of which we shall discuss in the 
next section. 

We now introduce some notation. We assume that the control- signal space U, the 
disturbance- signal space W, the error-signal space Z and the measurement signal 
space y consist of column vectors of given sizes nu, nyy, nz and ny, respectively, 
with entries from the quotient field K of A: 

We are given a plant G=[^ll %\l ] : W®U Z®y and seek to desi gn a controller 
K: y ^ U that stabilizes the system Yi{G,K) of Figured] as given in Section |2l 
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The various matrix entries Gij of G are now matrices with entries from K (rather 
than TZH°° as in the classical case) of compatible sizes (e.g., Gn has size nyy x nu) 
and if is a matrix over K of size nu x ny. Again vi and V2 are tap signals used to 
detect stability properties of the internal signals u and y. 

Just as was explained in Section [2] for the classical case, the system Yi{G,K) is 



well-posed if there is a well-defined map from 



to 



and this happens exactly 



when det(/ — G22K) ^ (where the determinant now is an element of A); when 



this is the case, the map from 



to 


" Z ' 

u 




. y . 



is given by 



e{G,K) 



where 0(G, K) is given by (|2.2p . We say that the system I](G, K) is internally stable 
if S(G, K) is well-posed and, in addition, if the map e(G, K) maps A"-^ qA"" A">' 
into A"^ © A"" © A"J', i.e., stable inputs w,vi,V2 are mapped to stable outputs 
z, u, y. Note that this is the same as the entries of S(G, K) being in A. 

To formulate the standard problem of i?°°-control, we assume that A is equipped 
with a positive-definite inner product making A at least a pre-Hilbert space with 
norm || • ||a; in the classical case, one takes this norm to be the i^-norm over the 
unit circle. Then we say that the system 'S(G,K) has performance if S(G, if) is 
internally stable and in addition the transfer function T^^, from w to z has induced 
operator norm bounded by some tolerance which we normalize to be equal to 1: 



sup{||z||^ 



klU-w < = 0, V2 = 0}< 1. 



We say that the system 'E,{G,K) has strict performance if in fact HT^^Hop < 1. 
The stabilization problem then is to describe all (if any exist) internally stabilizing 
controllers K for the given plant G, i.e., all K G IK""^"^ so that the associ- 
ated closed-loop system 'E,{G,K) is internally stable. The standard H°° -problem is 
to find all internally stabilizing controllers which in addition achieve performance 
||7ztu||op < 1. The strictly suboptimal iJ°°-problem is to describe all internally 
stabilizing controllers which achieve strict performance HT^u^Hop < 1. 

The ii°°-control problem for the special case where G22 = is the Model- 
Matching problem for this setup. With the same arguments as in Subsection l2.1l it 
follows that stabilizability forces Gn, G12 and G21 all to be stable (i.e., to have all 
matrix entries in A) and then K stabilizes exactly when also K is stable. 

3.1. Parametrization of stabilizing controllers in terms of a given stabiliz- 
ing controller. We return to the general case i.e.. 
Now suppose we have a stabilizing controller K e K""^"^. Set 



[/=(/- G22Kr^ and V = 

Then U G A"J'^"J', V G A""^"J', detC/ 7^ G 
Furthermore, 6(G, K) can then be written as 

Gil + Gi2t^G2i 
e{G,K) = e{G;U,V) := VG21 

UG21 



K{I 

L K -- 



G22K)-'. 
VU-^ and U 



(3.1) 

G22^^ = /. 



G12 + G12VG22 G12V 

I + VG22 V 
UG22 u 



(3.2) 



It is not hard to see that if [/ G A^v^'^y and V G A"""""*" are such that det [/ 7^ 0, 
U - G22V = I and 1321) is stable, i.e., in A("^+""+"y)^("^+""+"J'), then K = 

13 



VU ^ is a stabilizing controller. A dual result holds if we set 

U^{I-KG22y^ and V ^ {I - KG22y^ K. (3.3) 

In that case G A"" , i/ e A"" ^"J' , dot [7 7^ £ A, iiT = [/-i^, [/ - VG22 = I 
and we can write Q{G, K) as 



e{G,K) = Q{G]U,V) 



Gil + G12VG21 G12U G12V 

VG21 U _^ V ^ 

{I + G22V)G2l G22U I+G22V 



(3.4) 



while conversely, for any U G A"""""" and V £ A"" ''"J' withjiett/ ^ and 
U — VG22 — I and such that (13. 4|) is stable, we have that K — U^^V is a stabilizing 
controller. 

This leads to the following first-step more linear reformulation of the definition 
of internal stabilization. 

Theorem 3.1. A plant G defined by a transfer matrix G £ K^^z+^y)x(^w+nu) 
internally stabilizahle if and only if one of the following equivalent assertions holds: 

(1) There exists i = [^] e ^(.^u+ny)+ny ^-^^ detU j^O such that: 

(a) The block matrix p.2p is stable (i.e., has all matrix entries in A), and 

(b) [-G22 I]L = I. 

Then the controller K = VU^^ internally stabilizes the plant G and we 
have: 

U={I-G22K)-\ V = K{I -G22K)-\ 

(2) There exists L=[u -v]£ A""^(""+"y) with AetU ^ Q such that: 

(a) The block matrix (|3.4p is stable (i.e., has all matrix entries in A), and 

(b) LUJ:=[u-v][a{A-I- 

If this is the case, then the controller K = U internally stabilizes the 
plant G and we have: 

U = {I - KG22y\ V = {I-KG22)~^K. 

With this result in hand, we are able to get a parametrization for the set of all 
stabilizing controllers in terms of an assumed particular stabilizing controller. 

Theorem 3.2. 

(1) Let G K""^"y be a stabilizing controller for G G K("^+"y)^("w+""). 
Define — {I — G22K^:)~^ and = K{I — G22-f^*)^^- Then the set of 
all stabilizing controllers is given by 



K ={V,+Q){U, + G22Qy\ 
where Q G K""^"^ is an element of the set 



(3.5) 



G12 

I 
G22 



Q [G2I G22 l] G A.^^^^+^^+^y'l^^-^w+nu+ny) 



(3.6) 



such that in addition det(i7* + G22Q) 7^ 0. 

14 



(2) Let e K"«x«>> be a stabilizing controller for G G K("^+"i')^("w+""). 
Define ^ {I - K^G 22)^^ and = {I - K^G22y^K^. Then the set of 
all controllers is given by 

K =iU,+QG22y\K + Q), (3.7) 
where Q G is an element of the set 57 p.6p such that in addition 

det(f/, + QG22) 7^ 0. 

Moreover, if Q £ fl, that dct(C/, + G22Q) 7^ i/ and only z/det(C/, + QG22) 0, 
and the formulas (13. Sp and p.7p give rise to the same controller K . 



Proof. By Theorem 13. 1[ if ii' is a stabilizing controller for G, then K has the form 
K = VJJ-^ with L=[^]siS\n part (1) of Theorem [Q and then Q{G,K) is as in 
([3:21) . Similarly Q{G,K^) is given as e(G;t/*,V;) in ([321) with C/*,V; in place of 
[/, V. As by assumption 6(G; t/*, is stable, it follows that 0(G; [/, is stable 
if and only if e(G;C/,y)-e(G;;7*,K) is stable. Let Q = V-K; as t/ = I + G22V 
and = I + G22K, it follows that U = G22Q. From (|3.2p we then see that 
the stable quantity 9(G; U, V) — 6(G; C/,, K) is given by 



e(G;C/,y)-e(G;[/*,K) = 



G12 
/ 

G22 



g[G2i G22 /] 



Thus 

K = VU- 



(K + {V- + {U - t/,))-' = (K + g)((7* + G22Q)" 



where Q is an element of il such that det([/* + G22Q) 7^ 0. 

Conversely, suppose K has the form K = [V^ + Q){U^ 
and det([/» + G22Q) ^ 0. Define F = K + Q, t/ = C/, - 
checks that 



f G22(3)"^ where Q G 
G22Q- Then one easily 



e(G; [/,!/) = 9 (G;f/„V;) 



G12 
/ 

G22 



Q[G2i G22 /] 



is stable and 
[ -G22 / ] 



V 

U 



[ -G22 I ] 



[ -G22 I ] 



Q 
G22Q 



i+o = I. 



So = VU ^ stabilizes G by part (1) of Theorem l3.1l This completes the proof of 
the first statement of the theorem. The second part follows in a similar way by using 
the second statement in Theorem 13.11 and Q = V — V^,. Finally, since V = V and 
K = K, we find that indeed det{U^+G22Q) ^ if and only if det(?7* + QG22) ^ 0, 
and the formulas (|3.5p and (|3.7p give rise to the same controller K. □ 



The drawback of the parametrization of the stabilizing controllers in Theorem 
3.21 is that the set is not really a free-parameter set. By definition, Q G i7 if Q 
itself is stable (from the (1,3) entry in the defining matrix for the f2 in (|3.6p ). but, 
in addition, the eight additional transfer matrices 

G12QG21, G12QG22, G12Q, QG21, 

QG22, G22QG21, G22QG22, G22Q 

should all be stable as well. The next lemma shows how the parameter set Q, can 
in turn be parametrized by a free stable parameter A of size [nu + ny)x [uu + ny). 
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Lemma 3.3. Assume that G is stabilizable and that is a particular stabilizing 
controller for G. Let Q G K""^"^. Then the following are equivalent: 

(i) Q is an element of the set f2 in ()3.6p . 



(ii) 



Go 



Q [ G22 I ] is stable, 



(iii) Q has the form Q = LAL for a stable free-parameter A e A(""+">')^(""+"y), 
where L G A"" ^ and L e A(""+"J')^"i' are given by 



L=[{I- KM 



22 ) 



-{I -KM 



22 ) 



^K^ 



L 



-K4I - G22K,.)-^ 
{I-G22K,)-^ 



(3.8) 



Proof. The implication (i) (ii) is obvious. Suppose that A = [g^j] Q 1*^22 /] is 
stable. Note that 



LAL = [(/ - KM2 



= Q. 



-{L-KM22)-^K, 



22 



G22_ 

-K4L ~G22K,)-^ 
{I-G22K,)-^ 



Hence (ii) implies (iii). Finally assume Q ~ LAL for a stable A. To show that 
Q G r2, as A is stable, it suffices to show that 



Li := 



G12 

/ 
G22 



L is stable, and L2 := L [G21 G22 -f] is stable. 



Spelling out Li, using the definition of L from (|3.8p . gives 



G12 

/ 
G22 



[{I~KM22)-^ -{L-KM22)^^K,] 



We note that each of the six matrix entries of Li are stable, since they all occur 
among the matrix entries of 8(G, if*) (see (I2.2p ) and if* stabilizes G by assumption. 
Similarly, each of the six matrix entries of L2 given by 



L. 



-K4L-G22K, 



[G21 G22 i] 



(i - G22if*)^' 

is stable since if, stabilizes G. It therefore follows that Q G as wanted. 



□ 



We say that K stabilizes G22 if the map [Zl] 1-^ [y] in Figure [T] is stable, i.e., 
the usual stability holds with w = and z ignored. This amounts to the stability 
of the lower right 2x2 block in e(G, K): 

{L-KG22y^ {I-KG22)-^K 
G22{I - KG22)-^ I+G22{I-KG22)-^K 

The equivalence of (i) and (ii) in Lemma 13.31 implies the following result. 

Corollary 3.4. Assume that G is stabilizable. Then K stabilizes G if and only if 
K stabilizes G22- 
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Proof. Assume K^, G K""^"^ stabilizes G. Then in particular the lower left 2x2 
block in Q{G,K^) is stable. Thus stabilizes G22. Moreover, K stabilizes G22 if 
and only if K stabilizes G when we impose Gn = 0, G12 = and G21 = 0, that is, 
K is of the form ([33]) with and K as in Theorem El and Q e ]K""><"y is such 
that [1^22] Q I*^^^ ^] stable. But then it follows from the implication (ii) (i) 
in Lemma [3.31 that Q is in 17, and thus, by Theorem 13. 2[ K stabilizes G (without 
Gil = 0, G12 = 0, G21 = 0). □ 



Combining Lemma 13.31 with Theorem 13.21 leads to the following generalization 
of Theorem 12.11 giving a parametrization of stabilizing controllers without the as- 
sumption of any coprime factorization. 

Theorem 3.5. Assume that G G ]K("^+"y)^(""'+"") is stabilizable and that 
is one stabilizing controller for G. Define ~ {I — G22K^)~^ , V* — K^:{I — 
G22K^)-^, = {I - KM22)'^ and K = (/ - ii:*G22)"^^i:*. Then the set of all 
stabilizing controllers for G are given by 

K = {V, + Q)iU, + G22Q)"' = (C/* + QG22)"'(K + Q), 

where Q = LAL where L and L are given by 
of size {nu +ny) x [uk +ny) so that dct([/^ 
QG22) ^ 0. 



and A is a free stable parameter 
G22Q) =/= or equivalently det(?7* + 



3.2. The Youla-Kucera parametrization. There are two drawbacks to the para- 
metrization of the stabilizing controllers obtained in Theorem 13.51 namely, to find 
all stabilizing controllers one first has to find a particular stabilizing controller, and 
secondly, the map A i— s- Q given in Part (iii) of Lemma 13.31 is in general not one-to- 
one. We now show that, under the additional hypothesis that G22 admits a double 
coprime factorization, both issues can be remedied, and we are thereby led to the 
well known Youla-Kucera parametrization for the stabilizing controllers. 

Recall that G22 has a double coprime factorization in case there exist stable 
transfer matrices D, N, X, Y, D, N, X and Y so that the determinants of D, D, 
X and X are all nonzero (in A) and 



G2 



D-^iV = ND- 



1 


■ D 


-N' 




X 


N 




' Iny 


1 


-Y 


X 




Y 


b 




. /„„ 



(3.9) 



According to CoroUarv 13.41 it suffices to focus on describing the stabilizing con- 
trollers of G22 . Note that K stabilizes G22 means that 



(/-XG22)- 

G22(/-i^G22) 



{I-KG22)-^K 
I + G22{I - KG22)-^K 



is stable, or, by Theorem[321 that K is given by jSS]) or ^Ji) for some Q e ]K""^"J' 
so that [022] Q [ '^22 / ] is stable. 

In case G22 has a double coprime factorization Quadrat shows in [1201 Proposi- 
tion 4] that the equivalence of (ii) and (iii) in Lemma 13.31 has the following refine- 
ment. We provide a proof for completeness. 

Lemma 3.6. Suppose that G22 has a double coprime factorization (j3.9p . Let Q G 
j^«„xny_ y^jg^ ^^/j Q [G22 /] is stable if and only if Q = DAD for some A G 

^nuxny 
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Proof. Let Q = DAD for some A e A""^"^. Then 



/ 

G22 



Q [ G22 / ] 



QG22 Q 
G22QG22 G22Q 



DAN DAD 
NAN NAD 



Hence [ ] Q [ G22 i ] is stable. 

Conversely, assume that [g^^] Q[g22 i] is stable. Set A = D^^QD^^ . Then 
with X, Y, X and Y the transfer matrices from the coprime factorization (j3.9p we 
have 



A 



X -Y 



X -Y 



X -Y 



D 

N 


A[ N D ] 


■ -Y 
X 


DAN DAD 




■ -Y ' 


NAN NAD 




X 



QG22 Q 
G22QG22 G22Q 



~Y 
X 



Thus A is stable. 



□ 



Lemma 3.7. Assume that G22 admits a double coprime factorization p.9p . Then 
Kq is a stabilizing controller for G22 if o,nd only if there exist Xq G 
Yo e A""^"y, Xq £ A""^"" and Yq G A""><"y with det{Xo) + 0, det(Xo) ^ so 
that Kq = YqXq^ = X^'^Yq and 



D 


-N ' 




^0 


N 




■ Iny 


-Yo 


^0 . 




. ^0 


D 




. /„„ 



In particular, K = YX-^ = X-^Y is a stabilizing controller for G22, where 
X,Y,X,Y come from the double coprime factorization (j3.9p for G22- 



Proof. Note that if K is a stabilizing controller for G22, then, in particular, 

-1 



/ 

—G22 



-K 
I 



{I~KG22)-^ K{I-G22K)-^ 



(3.10) 



{I -G22K)-^G22 {I-G22K)-^ 

is stable. The above identity makes sense, irrespectively of K being a stabilizing 
controller, as long as the left hand side is invertible. Let X, Y, X and Y be the 
transfer matrices from the double coprime factorization. Set K = X^^Y = YX^^. 
Then we have 



X 


-Y ' 


-1 


' X 







-N 


D 







D 


^ ( 





/ 

-D-^N 





D-^ 



X 

-N 



-Y 
D 



-X-^Y 
I 



I 

—G22 



-K 
I 



Since X, D and 







-1 




X 


-Y 




D Y 


-N 


D 




_N X _ 



are stable, it follows that the right-hand 
side of (|3.10p is stable as well. We conclude that K = X^^Y = YX^^ stabilizes 

^22- 

Now let Kq be any stabilizing controller for G22- It follows from the first part 
of the proof that K — YX^^ — X^^Y is stabilizing for G22- Define V and U by 
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and V and U by 

a A e A""^"y so that 



Then, using Theorem l3.2l and Lemma l3.61 there exists 



Ka = {V + Q){U + = {U + QG22r\V + Q), 

where Q — DAD. We compute that 

(/ - G22K)-'^ = {I - D-^NYX-^)-^ = X{DX - NY)-^D = XD 

and 



(/ - KG22y^ = {I - X-^YND-^)-^ ^ D{XD - YN)-^X = DX. 



(3.11) 
(3.12) 



Thus 
Therefore 



V = YD, U = XD, V^DY, U = DX. 



and 



Set 



Ko 



(V + Q){U + = {YD + DKD){XD + NKD)-^ (3.13) 
{Y + bK){X + NK)-^ 

{U + QG22y\V + Q) = {DX + DANy\DY + DAD) (3.14) 
{X + AN)-\Y + AD). 



Yo^{Y + DA), Xo^{X + NA), Yo^{Y + AD), Xo = (X + AN). 
Then certainly det Xq ^ and dot Xq ^ 0, and we have 



D 


-N ' 




X 


N 




-% 


Xo . 






Yo D 






/ 


■ 




■ D 


-N' 




-A 


/ 




-Y 


X 



D ^-N 
-Y-AD X + AN 



X + NA N 
Y + DA D 



X N 
Y D 



I 
A / 



/ 
-A / 



/ 
A / 



/ 
/ 



□ 



Since any stabiUzing controller for G is also a stabilizing controller for G22, the 
following corollary is immediate. 

Corollary 3.8. Assume that G £ K("^+"y)^ is a stabilizable and that 
G22 admits a double coprime factorization. Then any stabilizing controller K of G 
admits a double coprime factorization. 

Lemma 3.9. Assume that G is stabilizable and that G22 admits a double coprime 
factorization. Then there exists a double coprime factorization (j3.9p for G22 so that 
DG21 and G12D are stable. 

Proof. Let ii' be a stabilizing controller for G. Then K is also a stabilizing con- 
troller for G22. Thus, according to Lemma 13.71 there exists a double coprime 
factorization ([33]) for G22 so that K = YX^^ = XY~^. Note that jSJ]) implies 
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that ^1 _Y ^ ^ I- In particular, DY = YD and NX = XN. Moreover, 
from the computations p. lip and (|3.12p we see that 

[I ~ G22K)-^ = XD and {I - KG22y^ = DX . 

Inserting these identities into the formula for 0(G, K), and using that K stabilizes 
G, we find that 

<d{G,K) 



Gil + G12YDG21 G12DX G12DY 
YDG21 DX DY^ 

XDG21 NX I + NY 



is stable. 



In particular 



G12DX G12DY 
G12DX G12DY 



is stable, and thus 



N^ 
-D 



= Gi2D{XN - YD) = G12D 



is stable. Similarly, since [ ]f 2g2i ] stable, we find that 



[^N D] 



YDG21 
XDG21 



{-NY + DX)DG2i = DG21 



is stable. 



□ 



We now present an alternative proof of Corollary 
a double coprime factorization. 



for the case that G22 admits 



Lemma 3.10. Assume that G is stabilizahle and G22 admits a double coprime 
factorization. Then K stabilizes G if and only if K stabilizes G22 ■ 

Proof. It was already noted that in case K stabilizes G, then K also stabilizes G22. 
Now assume that K stabihzes G22. Let Q e so that K is given by (|^ . 

It suffices to show that Q ^ fl, with defined by p.6p . Since G is stabilizahle, it 
follows from Lemma 13.91 that there exists a double coprime factorization (|3.9p of 
G22 so that DG21 and G12D are stable. According to Lemma [3.61 Q = DAD for 
some A e A""^"^. It follows that 



G12 

/ 
G22 



Q [ G21 G22 / ] 



G12D 
D ^ 
G22D 

G12D 
D 

N 



A [ DG21 DG22 D ] 



A [ DG21 N D ] 



is stable. Hence Q G fi. 



□ 



Combining the results from the Lemmas 13 . 61 13 . 71 and 13.10] with Theorem l3.2l and 
the computations (|3.13p and (|3.14p from the proof of Lemma 13.71 we obtain the 
Youla-Kucera paranietrization of all stabilizing controllers. 

Theorem 3.11. Assume that G € IK("^+"y)^(">^+"") is stabilizable and that G22 
admits a double coprime factorization l|3.9p . Then the set of all stabilizing con- 
trollers is given by 

K ={Y + DA){X + NA)-^ = (A + AN)-^{Y + AD), 
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where A is a free stable parameter from A""^"^^ such that det{X + iVA) ^ or 
equivalently det{X + AN) ^ 0. 

3.3. The standard _ff°°-problem reduced to model matching. We now con- 
sider the i/°°-problem for a plant G = [ ] : W ® ^ Z © J^, i.e., we seek a 
controller K : 3^ — > W so that not only ld{G,K) in (|2.2p is stable, but also 

IIGii + Gi2K{I - G22K)-'G2i\\op < 1. 

Assume that the plant G is stabilizable, and that X* : y U stabilizes G. 
Define C/* , , C/* and as in Theorem 13.21 We then know that all stabilizing 
controllers of G are given by 

K^{V, + Q){U, + G22Qy^ = [U, + QG22r\V, + Q) 

where Q G K""^"^ is any element of SI in p.6p . We can then express the transfer 
matrices U and V in (|3.ip in terms of Q as follows: 

U = {I - G22K)-' = {I - G22{V, - Q){U, - G22QrY' 
= {U, - G22Q){U, - G22Q ~ G22{V, - Q))-^ 
= -G22Q)(C/* 
- (C/*-G22Q), 
where we used that — G22^* — I, and 

V ^ KU ^V^- Q. 
Similar computations provide the formulas 

U = U, + QG22 and V = V, + Q 

for the transfer matrices U and V in p.3p . Now recall that Q{G,K) can be ex- 
pressed in terms of C/ and as in (|3.2p . It then follows that left upper block in 
0(G, i^T) is equal to 

Gn+Gi2K{I -G22K)-'^G2i = Gn + G12VG21 (3.15) 

= Gil + G12K.G21 — G12QG21. 

The fact that K^, stabilizes G implies that Gn := Gn + Gi2T4G2i is stable, and 
thus G12QG21 is stable as well. We are now close to a reformulation of the 
problem as a model matching problem. However, to really formulate it as a model 
matching problem, we need to apply the change of design parameter Q 1— > A defined 
in Lemma 13.31 or Lemma 13.61 in case G22 admits a double coprime factorization. 
The next two results extend the idea of Theorem l2.2l to this more general setting. 

Theorem 3.12. Assume that G £ K("^+">')^("^+"") is stabilizable and let K e 
K""^"y . Then K is a solution to the standard H°° problem for G if and only if 

K={V, + Q){U, + G22O)"' = (f/* + 0G22)"'(V; + 0) 

with Q = LAL, where L and L are defined by p.Sp . so that det(C/* + G22Q) 7^ 0, 
or equivalently det(C/* + QG22) 7^ 0, and A G Ai^u+ny)>i{nu+ny) ggi^tion to 

the model matching problem for Gn, G12 and G21 defined by 

Gil '■— Gil + G12KG21, G12 '■— G12L, G21 '■— LG21, 

i.e., so that 

||Gii + Gi2AG2i||op< 1. 
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Proof. The statement essentially follows from Theorem 13.51 and the computation 
p.lSp except that we need to verify that the functions Gn, 6*12 and G21 satisfy 
the conditions to be data for a model matching problem, that is, they should be 
stable. It was already observed that Gn is stable. The fact that G12 and G21 are 
stable was shown in the proof of Lemma 13.31 □ 

We have a similar result in case G22 admits a double coprime factorization. 

Theorem 3.13. Assume that G G k("^+">')^("w+"") is stabilizable and that G22 
admits a double coprime factorization p.9p . Let K G IK"^^"". Then K is a 
solution to the standard H°^' problem for G if and only if 

K={Y + DA){X + NA)-^ = iX + AAr)-i(y + AD), 

where A G A""''"^ so that det{X + NA) ^ 0, or equivalently dei{X + AN) ^ 0, is 
any solution to the model matching problem for Gn, G12 and G21 defined by 

Gil '■= Gil + Gi2i^£'G2i, G12 G12D, G21 := DG21, 

i.e., so that 

||Gn + Gi2AG2i||op< 1. 

Proof. The same arguments apply as in the proof of Theorem 13.121 except that in 
this case Lemma [3.91 should be used to show that G12 and G21 are stable. □ 

3.4. Notes. The development in Section [01 on the parametrization of stabilizing 
controllers without recourse to a double coprime factorization of G22 is based on 
the exposition of Quadrat [120] . It was already observed by Zames-Francis [140] 
that Q = K{I — G22K)~^ can be used as a free stable design parameter in case G22 
is itself already stable; in case G22 is not stable, Q is subject to some additional 
interpolation conditions. The results of |120j is an adaptation of this observation 
to the general ring-theoretic setup. The more theoretical papers [118t I122j give 
module-theoretic interpretations for the structure associated with internal stabiliz- 
ability. In particular, it comes out that every matrix transfer function G22 with 
entries in K has a double-coprime factorization if and only if A is a Bezout domain, 
i.e., every finitely generated ideal in A is principal; this recovers a result already 
appearing in the book of Vidyasagar [136j . A new result which came out of this 
module-theoretic interpretation was that internal stabilizability of a plant G22 is 
equivalent to the existence of a double-coprime factorization for G22 exactly when 
the ring A is projective-free, i.e., every submodule of a finitely generated free module 
over A must itself be free. This gives an explanation for the earlier result of Smith 
[130] that this phenomenon holds for the case where A is equal H°° over the unit 
disk or right-half plane. 

Earlier less complete results concerning parametrization of the set of stabilizing 
controllers without the assumption of a coprime factorization were obtained by 
Mori [102] and Sule |132j . Mori |103j also showed that the internal-stabilization 
problem can be reduced to model matching form for the general case where the 
plant has the full 2 x 2-block structure G = [%\\ %]^^ ] . 

Lemma 13.101 for the classical case is Theorem 2 on page 35 in [64] . The proof 
there relies in a careful analysis of signal-fiow diagrams; we believe that our proof 
is more direct. 
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4. Feedback control for linear time-invariant multidimensional 

SYSTEMS 

4.1. Multivariable frequency-domain formulation. The most obvious multi- 
variable analogue of the classical single-variable setting considered in the book of 
Francis is as follows. We take the underlying field to be the complex num- 
bers C; in the engineering applications, one usually requires that the underlying 
field be the reals R, but this can often be incorporated at the end by using the 
characterization of real rational functions as being those complex rational func- 
tions which are invariant under the conjugation operator s{z) >—>■ s{z). We let 
D'^ = {z = (zi, . . . , Zd) ■ \zk\ < 1} be the unit polydisk in the c?-dimensional com- 
plex space and we take our ring A of stable plants to be the ring C{z)s of all 
rational functions s{z) = in d variables (thus, p and q are polynomials in the d 
variables zi, . . . ,Zd where we set z — (zi, . . . , z^)) such that s(z) is bounded on the 
polydisk D''. The ring C[z] of polynomials in d variables is a unique factorization 
domain so we may assume that p and q have no common factor (i.e., that p and 
q are relatively coprime) in the fractional representation s = ^ for any element of 
C(zi, . . . , zj). Unlike in the single- variable case, for the case c? > 1 it can happen 
that p and q have common zeros in even when they are coprime in C[z] (see 
[138] for an early analysis of the resulting distinct notions of coprimeness). It turns 
out that for d > 3, the ring C(z)s is difficult to work with since the denominator q 
for a stable ring element depends in a tricky way on the numerator p: if s g C{z)s 
has coprime fractional representation s — ^, while it is the case that necessarily 
q has no zeros in the open polydisk D'', it can happen that the zero variety of q 
touches the boundary as long as the zero variety of p also touches the same 
points on the boundary in such a way that the quotient s — ^ remains bounded on 

D''. Note that at such a boundary point (, the quotient s = p/q has no well-defined 
value. In the engineering literature (see e.g. [45 ) 11311 184]). such a point is known 
as a nonessential singularity of the second kind. 

To avoid this difficulty, Lin [92l [93] introduced the ring C(z)ss of structured 
stable rational functions, i.e., rational functions s G C(z) so that the denominator 
q in any coprime fractional representation s — ^ for s has no zeros in the closed 

polydisk ID)'^. According to the result of Kharitonov-Torres-Muhoz [84], whenever 
s = 2 £ C(z)s is stable in the ffist (non-structured) sense, an arbitrarily small 
perturbation of the coefficients of q may lead to the perturbed q having zeros in 
the open polydisk B'' resulting in the perturbed version s — ^ oi s being unstable; 
this phenomenon does does not occur for s E C{z)ss, and thus structured stable 
can be viewed just as a robust version of stable (in the unstructured sense). Hence 
one can argue that structured stability is the more desirable property from an 
engineering perspective. In the application to delay systems using the systems- 
over-rings approach [46l [85l [83], on the other hand, it is the collection C{z)ss of 
structurally stable rational functions which comes up in the first place. 

As the ring A = C{z)ss is a commutative integral domain, we can apply the 
results of Section [3] to this particular setting. It was proved in connection with 
work on systems-over-rings rather than multidimensional systems (see |461 183] ) 
that the ring C{z)ss is projective-free. As pointed out in the notes of Section [3] 
above, it follows that stabilizability of G22 is equivalent to the existence of a double 
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coprime factorization for the plant 6*22 (see [119) ). thereby setthng a conjecture of 
Lin [521[S31[M]. We summarize these results as follows. 

Theorem 4.1. Suppose that we are given a system G — [q^^ gH] over the quo- 
tient field Q{<C{z)ss) of the ring C{z)ss of structurally stable rational functions in 
d variables. If there exists a controller K = YX~^ = X~^Y which internally sta- 
bilizes G, then G22 has a double coprime factorization and all internally stabilizing 
controllers K for G are given by the Youla-Kucera parametrization. 

Following Subsection l3.3[ the Youla-Kucera parametrization can then be used to 
rewrite the iJ °°-problem in the form of a model-matching problem: Given Ti , T2 , Ts 
equal to matrices over <C{z)ss of respective sizes nz x nyv, 't-w x and ny x nw, 
find a matrix A over <C{z)ss of size ny x ny so that the affine expression S given 
by 

S^Ti+T^An (4.1) 

has supremum norm at most 1, i.e., ||S'||oo = max{||S'(z)|| : z G B'*} < 1. 

For mathematical convenience we shall now widen the class of admissible so- 
lutions and allow Ai, . . . , Aj to be in the algebra i7°°(I!)'') of bounded analytic 
functions on D''. The unit ball of iJ°°(D'*) is the set of all holomorphic functions S 
mapping the polydisk 3'^ into the closed single- variable unit disk D C C; we denote 
this space by Sd, the d-variable Schur class. While Ti, T2 and T3 are assumed to 
be in C{z)ss, we allow A in (gH]) to be in H°°(p''). 

Just as in the classical one- variable case, it is possible to give the model-matching 
form l|4.ip an interpolation interpretation, at least for special cases (see [73l [74l [32] ) . 
One such case is where nyy = nz = ny ~ 1 while nu — J . Then T\ and T3 are 
scalar while T2 = [Tb.i •■■ T2,.7 ] is a row. Assume in addition that T3 = 1. Then the 
model-matching form (|4.ip collapses to 

5 - Ti + T21A1 + ■ • • + T2,7Aj (4.2) 

where Ai,...Aj are J free stable scalar functions. Under the assumption that 

the intersection of the zero varieties of T^.x, . . . , T2,j within the closed polydisk D 
consists of finitely many (say iV) points 

and if we let wi, . . . , wat be the values of T\ at these points 

wx = Ty{z{), . . . , WW = ri(z7v), 

then it is not hard to see that a function S G C(z)ss has the form (|4.2p if and only 
if it satisfies the interpolation conditions 

S'(z,) =Wi fori = l,...,iV. (4.3) 

In this case the model-matching problem thus becomes the following finite-point 
Nevanlinna-Pick interpolation problem over B'^: find S G C{z)ss subject to |>5'(z)| < 
1 for all z G B'' which satisfies the interpolation conditions (|4.3p . Then the d- 
variable i/°°-Model-Matching problem becomes: find S E Sd so that S{zi) = wi 
fori = l,...,N. 

A second case (see [32]) where the polydisk Model-Matching Problem can be 
reduced to an interpolation problem is the case where T2 and T3 are square (so 
nz = nu and ny = nw) with invertible values on the distinguished boundary of 
the polydisk; under these assumptions it is shown in [32] (see Theorem 3.5 there) 
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how the model-matching problem is equivalent to a bitangential Nevanlinna-Pick 
interpolation problem along a subvariety, i.e., bitangential interpolation conditions 
are specified along all points of a codimension-1 subvariety of D'' (namely, the union 
of the zero sets of detT2 and detTa intersected with D*^). For d = 1, codimension-1 
subvarieties are isolated points in the unit disk; thus the codimension-1 interpola- 
tion problem is a direct generalization of the bitangential Nevanlinna-Pick inter- 
polation problem studied in [26l [58l [62] . However for the case where the number 
of variables d is at least 3, there is no theory with results parallel to those of the 
classical case. 

Nevertheless, if we change the problem somewhat there is a theory parallel to 
the classical case. To formulate this adjustment, we define the d-variable Schur- 
Agler class SAd to consist of those functions S analytic on the polydisk for which 
the operator S{Xi, . . . , X^) has norm at most 1 for any collection Xi, . . . , X^ of 
d commuting strict contraction operators on a separable Hilbert space /C; here 
S{Xi, . . . , Xd) can be defined via the formal power series for S: 

S{X,,...,Xd)^ SnX'\ if5(z)= 

where we use the standard multivariable notation 

n = (ni, ...,nd)e 1%, X" = X'l' ■ ■ ■ Xy and z" = zj^^ • • • z^^ 

For the cases d = 1, 2, it turns out, as a consequence of the von Neumann inequality 
or the Sz.-Nagy dilation theorem for d = 1 and of the Ando dilation theorem [T7] 
for d — 2 (see |109[ [34] for a full discussion) , that the Schur- Agler class SAd and 
the Schur class Sd coincide, while, due to an explicit example of Varopoulos, the 
inclusion SAd C Sd is strict for d > 3. 

There is a result due originally to Agler [1] and developed and refined in a 
number of directions since (see [31 135j and jjj for an overview) which parallels the 
one-variable case; for the case of a simple set of interpolation conditions (|4.3p the 
result is as follows: there exists a function S in the Schur-Agler class SAd which 
satisfies the set of interpolation conditions S{zi) = wi for i = I, . . . , N if and only 
if there exist d positive semidefinite matrices . . . , P'-'') of size N x N so that 

d 

fc=l 

For the case d — I, the Pick matrix P = [ \^^]^- ] ^ is the unique solution of 
this equation, and we recover the classical criterion P > for the existence of 
solutions to the Nevanlinna-Pick problem. There is a later realization result of 
Agler |2J (see also O [35]): a given holomorphic function S is in the Schur-Agler 
class SAd{^il^,y)) if and only if S has a contractive Givone-Roesser realization: 
S{z) = D + C{I - Z{z)A)-^Z{z)B where [^g] : {^f^^Xk ® U) ^ {(Bf^^Xk (B y) 

is contractive with Z(z) = 

L ^dlXj J 

Direct application of the Agler result to the bitangential Nevanlinna-Pick inter- 
polation problem along a subvariety, however, gives a solution criterion involving an 
infinite Linear Matrix Inequality (where the unknown matrices have infinitely many 
rows and columns indexed by the points of the interpolation-node subvariety) — see 
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[321 Theorem 4.1]. Alternatively, one can use the polydisk Commutant Lifting The- 
orem from |31] to get a solution criterion involving a Linear Operator Inequality 
[521 Theorem 5.2]. Without further massaging, either approach is computationally 
unattractive; this is in contrast with the state-space approach discussed below. In 
that setting there exists computable sufficient conditions, in terms of a pair of LMIs 
and a coupling condition, that in general are only sufficient, unless one works with 
a more conservative notion of stability and performance. 

4.2. Multidimensional state-space formulation. The starting point in this 
subsection is a quadruple {A, B, C, D} consisting of operators A, B, C and D 
so that [ ^ B ] ■ [ w%u ] ~^ [ z®y ] ^'^d a partitioning X = Xi® ■ ■ ■ ® Xd oi the space 
X. Associate with such a quadruple {A^ B, C, D} is a linear state-space system S 
of Givone-Roesser type (see [57] ) that evolves over and is given by the system 
of equations 



" xi{n+ei) ' 




' xi{n) ' 




= A 








.Xd(n). 



Bu{n) 



in e z^). 



(4.4) 



y{n) = Cx{n) + Du{n) 
with initial conditions a specification of the state values XkC^j-n^tjCj) for t 



{ti, . . . ,td) e Z"^ subject to tk — where k — 1 

' xi(n) 

unit vector in and x{n) — ; 

-Xd(n) 

operator. Moreover, the block operator matrix [ ^ |^ ] is referred to as the syst 



, d. Here Cfe stands for the k-ih 
We call X the state-space and A the state 



matrix. 

Following [81] , the Givone-Roesser system (|4.4p is said to be asymptotically stable 
in case, for zero input u{n) — for n £ and initial conditions with the property 

d 

sup ||a^/c(^^ < oo for = 1, . . . , d, 

the state sequence x satisfies 

sup ||a:(n)|| < oo and lim ||2:(n)|| ~ 0, 



where n — > oo is to be interpreted as min{ni, . . . nd\ — > oo when n = (ni, . . . , Ud) G 
Z^. 

With the Givone-Roesser system (|4.4[) we associate the transfer function G{z) 
given by 

G{z) = D + C{I - Z{z)A)-^Z{z)B, (4.5) 
defined al least for z G with ||z|| sufficiently small, where 



Z{z) 



ZdlXd 



(zgC^). 



(4.6) 



We then say that {A, i?,C, Z?} is a (state-space) realization for the function G, or 
if G is not specified, just refer to {A, B, G, D} as a realization. The realization 
{A, B, C, D}, or just the operator A, is said to be Hautus-stable in case the pencil 
/ — Z{z)A is invertible on the closed polydisk D . 
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Here we only consider the case that X is finite dimensional; then the entries 
of the transfer function G are in the quotient field Q{<C{z)ss) of C(z)ss and are 
analytic at 0, and it is straightforward to see that G is structurally stable in case 
G admits a Hautus-stable realization. For the case d = 2 it has been asserted in 
the literature [STJ Theorem 4.8] that asymptotic stability and Hautus stability are 
equivalent; presumably this assertion continues to hold for general d> \ but we do 
not go into details here. 

Given a realization {A, B, C, D} where the decomposition X ~ Xi Q) ■ ■ ■ (B is 
understood, our main interest will be in Hautus-stability; hence we shall say simply 
that A is stable rather than Hautus-stable. 

As before we consider controllers K in Q{C{z) 
assume to be given by a state-space realization 



of size ny x uk that we also 

K{z) = Dk + Ck{I~ Zk{z)Ak)-'Zk{z)Bk (4.7) 

with system matrix [ ] : ['^y] — > [ '^f ] i a- decomposition of the state-space 
Xk = Xi^K © ■ • ■ © Xd.K and Zk{z) defined analogous to Z{z) but with respect to 
the decomposition of Xk- We now further specify the matrices B, C and D from 
the realization {A, B, C, D} as 



B=[ Bi B2] 



C 



Ci 
C2 



D 



D21 



D12 
D22 



compatible with the decompositions Z 



U. 



(4.8) 

We can then form the 



3^ and W 

closed loop system Gd = S(G, i^T) of the two transfer functions. The closed loop 
system Gd = S(G, K) corresponds to the feedback connection 



A 
C2 

subject to 



D21 





X 




X 




w 




z 




u 







' Ak 


Bk ' 




XK 




XK 




Dk . 




UK 




VK 



B2 
D12 
D22 



Z{z)x, xk = Zk{z)xk, UK = y and vk = 



This feedback loop is well-posed exactly when I — D22D k is invertible. Since, under 
the assumption of well posedness, one can always arrange via a change of variables 
that D22 = (cf., [78]), we shall assume that D22 = for the remainder of the 
paper. In that case well-posedness is automatic and the closed loop system Gd 
admits a state-space realization 



Gd{z) 
with system matrix 



Dd + C'dil - Zdiz)Ad)-'Zdiz)Cd 



Ad 
Cd 



Bd 
Dd 



A- 


\- B2DKC2 


B2CK 


B1 + B2DKD21 




BkC2 


Ak 


BkD21 


. ^1- 


¥D^2DkC2 


DuCk 


Dn + D12DKD21 _ 



(4.9) 



(4.10) 



and 



Zd{z) 



[z e C^). 



Z{z) 
Zk{z) 

The state-space (internal) stabilizability problem then is: Given the realization 
{A, B,C, D} find a compatible controller K with realization {Ak, Bk,Ck, Dk} 
so that the closed-loop realization {Ad, Bd,Cd, Dd} is stable, i.e., so that / — 
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Zci{z)Aci is invertible on the closed polydisk D . We also consider the strict state- 
space -problem: Given the realization {A, S, C, 13}, find a compatible con- 
troller K with realization {AktBktCktDk} so that the closed loop realization 
{Aci, BcuCci, Del} is stable and the closed-loop system Gd satisfies ||G'c;(z)|| < 1 
for all z G B'^. 



State-space stabilizability. In the fractional representation setting of Section[3]it 
took quite some effort to derive the result: "If G is stabilizable, then K stabilizes G 
if and only \i K stabilizes G22" (see Corollary 13.41 and Lemma FB.lOp . For the state- 
space stabilizability problem this result is obvious, and what is more, one can drop 
the assumption that G needs to be stabilizable. Indeed, G22 admits the realization 
{A, B2, C2, 0} (assuming D22 = 0), so that the closed-loop reahzation for I]{G22,K) 
is equal to {Aci, 0,0,0}. In particular, both closed-loop realizations have the same 
state operator Ad, and thus K with reahzation {Ak,Bk,Ck,Dk} stabihzes G22 
if and only if K stabilizes G, without any assumption on the stabilizability of G. 

The state-space stabilizability problem does not have a clean solution; To discuss 
the partial results which exist, we first introduce some terminology. 

Let {A, B, C, D} be a given realization as above with decomposition of B, C 
and D as in (|4.8|) . The Givone-Roesser output pair {C2, A} is said to be Hautus- 



detectable if the block-column matrix 



-Z{z)A 
C2 



is of maximal rank nx (i-C, is left 



invertible) for all z in the closed polydisk D''. We say that {(72, A} is operator- 
detectable in case there exists an output-injection operator L : ^ — > <Y so that 
A -|- LC2 is stable. Dually, the Givone-Roesser input pair {A, B2} is called Hautus- 
stabilizable if it is the case that the block-row matrix [/ — AZ{z) B2] has maximal 

rank nx (i.e., is right invertible) for all z £ and operator- stabilizable if there 
is a state-feedback operator F: X U no that A -\- B2F is stable. Notice that 
both Hautus-detectability and operator-detectability for the pair (G, A) reduce to 
stability of A in case G = 0. A similar remark applies to stabilizability for an input 
pair {A,B). 

We will introduce yet another notion of detectability and stabilizability shortly, 
but in order to do this we need a stronger notion of stability. We first define V to 
be the set 

P=| •. : X,:X,^ X„ i = l,,...,d^, (4.11) 

which is also equal to the commutant of {Z{z) : z G Z''} in the G*-algebra of 
bounded operators on X. We then say that the realization {A, B, G, D}, or just A, is 
scaled stable in case there exists an invertible operator Q ^ V so that ||(5~^A(5|| < 1, 
or, equivalently, if there exists a positive definite operator X (notation X > 0) in I? 
so that AX A* — X < 0. To see that the two definitions coincide, take either X = 
QQ* G or, when starting with X > 0, factor X as X = QQ* for some Q G 13. It 
is not hard to see that scaled stability implies stability. Indeed, assume there exists 
an invertible Q G 23 so that \\Q-^AQ\\ < 1. Then Z{z)Q-'^AQ = Q-^Z{z)AQ is a 

strict contraction for each z G D , and thus Q^^{I — Z{z)A)Q = I — Z{z)Q^^AQ 
is invertible on D . But then / — Z{z)A is invertible on ED as well, and A is stable. 
The converse direction, even though asserted in |11H [95] , turns out not to be true 
in general, as shown in [16] via a concrete example. The output pair {G2, A} is then 
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said to be scaled- detectable if there exists an output-injection operator L : 3^ — > <Y so 
that A + LC'2 is scaled stable, and the input pair {A, B2} is called scaled- stahilizable 
if there exists a state- feedback operator F : X so that A-f i?2-F is scaled stable. 

While a classical result for the 1-D case states that operator, Hautus and scaled 
detect ability, as well as operator, Hautus and scaled stabilizability, are equivalent, 
in the multidimensional setting considered here only one direction is clear. 

Proposition 4.2. Let {A, B be a given realization as above with decompo- 
sition of B, C and D as in (|4.8p . 

(1) If the output pair {C 2 T A} is scaled-detectable, then {C2,A} is also operator- 
detectable. If the output pair {C2, A} is operator- detectable, then {C2,A} 
is also Hautus- detectable. 

(2) If the input pair {A, B2} is scaled-stabilizable, then {A, B2} is also operator- 
stabilizable. If the input pair {A, B 2} is operator- stabilizable, then {A, B2} 
is also Hautus- stabilizable. 

Proof. Since scaled stability is a stronger notion than stability, the first implications 
of both (1) and (2) are obvious. Suppose that L: y ^ X is such that A -\- LC2 is 
stable. Then 



[/ -Z{z)L] 



- z{z)A 
C2 



= 1 - Z{z){A + LC2) 



is invertible for all z G E' ' from which it follows that {C2, A} is Hautus-detectable. 
The last assertion concerning stabilizability follows in a similar way by making use 
of the identity 

/ 

[-FZ{z)_ 

□ 



[l-AZ{z) B2] 



I-{A + B2F)Z{z). 



The combination of operator-detectability together with operator-stabilizability 
is strong enough for stabilizability of the realization {A, B, C, D} and we have the 
following weak analogue of Theorem [ 



Theorem 4.3. Let {A, B, C, D} be a given realization as above with decomposition 
of B, C and D as in ()4.8p (with D22 — 0). Assume that {C2,A} is operator- 
detectable and {A, B2} is operator- stabilizable. Then {A, B, C, D} is stabilizable. 
Moreover, in this case one stabilizing controller is K ~ {Ak, Bk,Ck, Dk} where 

\A + B2F + LC2 -L-\ 



'Ak 


Bk 




Ck 


Dk_ 





B2F 
F 







where L : y X and F : X 
A + FB2 are stable. 



U are any operators chosen such that A + LC2 ana 



Proof. It is possible to motivate these formulas with some observability theory (see 
|57| ) but, once one has the formulas, it is a simple direct check that 



'A^i 

Ccl 



Bel 
Del 



A + B2DKC2 

BkC2 

A 

-LC2 



B2CK 
Ak 



A 



B2F 
B2F + LC2 
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It is now a straightforward exercise to check that this last matrix can be put in the 
triangular form ["^^^ a+^2f] ^'^^ ^ sequence of block-row/block-column similar- 
ity transformations, from which we conclude that Ad is stable as required. □ 

Remark 4.4. A result for the systems-over-rings setting that is analogous to that 
of Theorem 14.31 is given in [5S] . There the result is given in terms of a Hautus- 
type stabilizable/detectable condition; in the systems-over-rings setting, Hautus- 
detectability/stabilizability is equivalent to operator-detectability/stabilizability (see 
Theorem 3.2 in [83]) rather than merely sufficient as in the present setting (see 
Proposition 14.21 above). The Hautus-type notions of detectability and stabilizabil- 
ity in principle are checkable using methods from [80] : see the discussion in [83l page 
161]. The weakness of Theorem 14. 31 for our multidimensional setting is that there 
are no checkable criteria for when {C2, A} and {A, B2} are operator-detectable and 
operator-stabilizable since the Hautus test is in general only necessary. 

An additional weakness of Theorem 14.31 is that it goes in only one direction: 
we do not assert that operator-detectability of {C2,A} and operator-stabilizability 
for {A, B2} is necessary for stabilizability of {A, B, C, D}. These weaknesses prob- 
ably explain why apparently this result does not appear explicitly in the control 
literature. 

While there are no tractable necessary and sufficient conditions for solving the 
state-space stabilizability problem available, the situation turns out quite differently 
when working with the more conservative notion of scaled stability. The following 
is a more complete analogue of Theorem 12.41 combined with Theorem 12.31 

Theorem 4.5. Let {A, B, C, D} be a given realization. Then {A, B, C, D} is scaled- 
stabilizable, i.e., there exists a controller K with realization {Ak, Bk,Ck, Dk} so 
that the closed loop state operator Ad is scaled stable, if and only if the input pair 
{A, B2} is scaled operator-stabilizable and the output pair {C2, ^4} is scaled operator- 
detectable, i.e., there exist matrices F and L so that A-\-B2F and A-\-LC2 are scaled 
stable. In this case the controller K given by (|4.12|) solves the scaled- stabilization 
problem for {A, B, C, D}. Moreover: 

(1) The following conditions concerning the input pair are eguivalent: 

(a) {A, B2} is scaled operator-stabilizable. 

(b) There exists Y E V satisfying the LMIs: 

B2,^{AYA* - Y)Bl^ < 0, y > (4.13) 

where B2.i_ any infective operator with range equal to Keri?2- 

(c) There exists Y E V satisfying the LMIs 

AY A* - Y - B2B; < 0, y > 0. (4.14) 

(2) The following conditions concerning the output pair are equivalent: 

(a) {C2A, } is scaled operator-detectable. 

(b) There exists X E T> satisfying the LMIs: 

Cl^{A*XA - X)C2,± < 0, a: > 0. (4.15) 

where C2.j_ any infective operator with range equal to KerC2. 

(c) There exists X eT> satisfying the LMIs 



A*XA - A - CIC2 < 0, A > 0. 
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(4.16) 



One of the results we shall use in the proof of Theorem 14.51 is known as Finsler's 
lemma [61 , which also plays a key role in [HHllZSj- This result can be interpreted as 
a refinement of the Douglas lemma [51j which is well known in the operator theory 
community. 

Lemma 4.6 (Finsler's lemma). Assume R and H are given matrices of appropriate 
size with H = H* . Then there exists a ^ > so that fiR*R > H if and only if 
R*^HR± < where R± is any injective operator with range equal to keri?. 

Finsler's lemma can be seen as a special case of another important result, which 
we shall refer to as Finsler's lemma II. This is one of the main underlying tools in 
the proof of the solution to the °°-problem obtained in [Ml HH] • 

Lemma 4.7 (Finsler's lemma II). Given matrices R, S and H of appropriate sizes 
with H = H* , the following are equivalent: 



(i) There exists a matrix J so that _ff + [ i?* S** ] 



with ranges equal to kcr R and ker S , respectively. 



' 


J* ' 




■ R ' 


J 







S 


are 


injec 


tive 



The proof of Finsler's Lemma II given in [66l uses only basic linear algebra and 
is based on a careful administration of the kernels and ranges from the various 
matrices. In particular, the matrices J in statement (i) can actually be constructed 
from the data. We show here how Finsler's lemma follows from the extended 



Proof of lemma \4.6\ using Lemma \4-. T\ Apply Lemma 14.71 with R = S. Then (ii) 
reduces to R*^HR± < 0, which is equivalent to the existence of a matrix J so that 
K ~ -(J* + J) satisfies R*KR > H. Since for such a matrix K we have K* — K, 
it follows that R*KR > H holds for K = ^il as long as ^il > K. □ 

With these results in hand we can proof Theorem 14.51 

Proof of Theorem \4-5\ We shall first prove that scaled stabilizability of {A, B, C, D} 
is equivalent to the existence of solutions X and Y in V for the LMIs (I4.15|) and 
(|4.13|) . Note that Ad can be written in the following affine way: 



Art = 



' A 


■ 




■ 


B2 ' 




■ Ak 


Bk ' 







I ' 








+ 


/ 









Dk . 




C2 






(4.17) 



Now let Xci : X Xk be an invertible matrix in Vd, where Vd stands for the 
commutant of {Zd{z): z G Z'^}. Let X be the compression of Xd to X and Y 
the compression of X^^^ to 
particular, X > and Y > 0. Then AdXdAd 

-x: ^ 



Then X,Y e V. Assume that Xd > 0. Thus, 
Xd < if and only if 



A* 



Ad 
—Xd 



< 0. 



(4.18) 



Now define 



H ^ 



A* 
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—Xd. 
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, R* = 









, s* = 


B2 
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J 



and 

■ Ak Bk 
Ck Dk 

Note that iJ, R and S are determined by the problem data, while J amounts to 
the system matrix of the controller to be designed. Then 



A*, 



Act 



H+[ R* S* ] 
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J* ■ 




■ R ' 


J 







s 



(4.19) 

Thus, by Finsler's lemma II, the inequality (I4.18P holds for some J = [c^ dZI 
and only if R*^HR± < and S'^HS± < 0, where without loss of generality we can 
take 



Ra 



I 











/ 























and Si 












B2,± 











/ 











/ 



with C2,± and -B2,_l as described in part (b) of statements 1 and 2. Writing out 
RIHR^ we find that R\HR^ < if and only if 



[ Cl^A* ] 



AC2,i_ 


-C2^XC2. 



< 



which, after taking a Schur complement, turns out to be equivalent to 

AC2,± 




Cl^{A*XA - X)C2.± = [ Cl^A* ] X, 



C;^XC2,± < 0. 



A similar computation shows that S*^HS± < is equivalent to B2^±{AYA* — 
Y)B2 < 0. This proves the first part of our claim. 



For the converse direction assume we have X and Y inV satisfying (|4.15p - (|4.13p . 
Most of the implications in the above argumentation go both ways, and it suffices to 
prove that there exists an operator Xd on X(BXk in 'DcU with Xk an arbitrary finite 
dimensional Hilbert space with some partitioning Xk = Xk,i © ■ ■ • ® Xx^d, so that 
Xci > and X and Y are the compressions to X oi Xd and X~i ^, respectively. Since 
()4.15p - ()4.13p hold with X and Y replaced by pX and pY for any positive number 
p, we may assume without loss of generality that ['/ y] > 0. The existence of the 
required matrix X^i can then be derived from Lemma 7.9 in [57] (with hk = n). 
To enforce the fact that Xci be in Dci we decompose X = diag(Xi, . . . ,Xd) and 
Y — diag(Yi, . . . , Yd) as in (|4.1ip and complete Xi and Yi to positive definite 
matrices so that [ * ] ^ = [ ^* * ] • 

To complete the proof it remains to show the equivalences of parts (a), (b) 
and (c) in both statements 1 and 2. The equivalences of the parts (b) and (c) 
follows immediately from Finsler's lemma with R = B2 (respectively, i? = C|) and 
H = AY A* - Y (respectively, H = A*XA - X), again using that X in (|4T5| can 
be replaced with pX (respectively, Y in (|4.13p can be replaced with pY) for any 
positive number p. 

We next show that (a) is equivalent to (b) for statement 1; for statement 2 the 
result follows with similar arguments. Let F : X ^ U, and let X G 1) be positive 
definite. Taking a Schur complement it follows that 



{A* + F*B;)X{A + B2F) -X <Q 
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(4.20) 



if and only if 



A 



'+F*BX 



A + B2F 

-X 



< 0. 



Now write 
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A + B2F 

-X 
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■ B2 
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F ' 




■ B* 
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-X 
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/ 




p* 
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Thus, applying Finsler's lemma II with 

^X-i A 
A* -X 



H = 



S=[0 I ] and J = i^, (4.21) 



, R^[BI ] 

we find that there exists an F so that (|4.20p holds if and only if 

RlHRi_ < and SlHS± < 

with now i?^ = [^p"^ ^] and S± ~ [q]. The latter inequality is the same as 
—X~^ < and thus vacuous. The first inequality, after writing out R*^HR^, turns 
out to be 

A*B2.i_ ~X 



<0, 



which, after another Schur complement, is equivalent to B2 i_{AX ^A* — X ^)i?2,i- 

'□ 

Since scaled stability implies stability, it is clear that finding operators F and L 
wit A + B2F and A + LC2 scaled-stable implies that A + B2F and A + LC2 are also 
stable. In particular, having such operators F and L we find the coprime factor- 
ization of G22 via the functions in Theorem 14. 31 While there are no known tractable 
necessary and sufficient conditions for operator-detectability/stabilizability, the LMI 
criteria in parts (iii) and (iv) of Theorem 14.51 for the scaled versions are considered 
computationally tractable. Moreover, an inspection of the last part of the proof 
shows how operators F and L so that A + B2F and A + LC2 are scaled stable can be 
constructed from the solutions X and Y from the LMIs in (|4.13p - (|4.16p : Assume 
we have X,Y e V satisfying (|4T3)) - ((4T6l) . Define H, R and S as in (g^, and 
determine a J so that -ff + [^J* S* ] [ j ] [f] < 0; this is possible as the proof of 
Finsler's lemma II is essentially constructive. Then take F = J. In a similar way 
one can construct L using the LMI solution Y. 

Stability versus scaled stability, /i versus /2. We observed above that the 
notion of scaled stability is stronger, and more conservative than the more intuitive 
notions of stability in the Hautus or asymptotic sense. This remains true in a more 
general setting that has proved useful in the study of robust control [SSI [571 1107] 
and that we will encounter later in the paper. 

Let A be a bounded linear operator on a Hilbert space X. Assume that in 
addition we are given a unital C*-algebra A which is realized concretely as a subal- 
gebra of C{X), the space of bounded linear operators on X. The complex structured 
singular value /iA(^) of A (with respect to the structure A) is defined as 

^^^^^ " inf{CT(A) : A e A, I - AA is not invertible} ' ^"^'^^^ 
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Here a{M) stands for the largest singular value of the operator M. Note that this 
contains two standard measures for A: the operator norm ||j4|| if we take A — C{JY), 
and p{A), the spectral radius of A, if we take A ~ {XIx ■ A S C}; it is not hard to 
see that for any unital C*-algebra A we have p{A) < fMA{A) < \\A\\. See [107] for 
a tutorial introduction on the complex structured singular value and [60j for the 
generalization to algebras of operators on infinite dimensional spaces. 

The C*-algebra that comes up in the context of stability for the iV-D systems 
studied in this section is A = {Z{z) : z G C^}. Indeed, note that for this choice of 
A we have that A is stable if and only if ^^{A) < 1. 

In order to introduce the more conservative measure for A in this context, we 
write 2?A for the commutant of the C*-algebra A in C{X). We then define 

Ma(^) — inf{7: < 7 for some invertible Q G 2?a} 

= inf{7: AX A* - 7X < for some X e Pa, ^ > 0}. (4.23) 

The equivalence of the two definitions again goes through the relation between X 
and Q via X = Q*Q. It is immediate that with A = {Z{z) : z G C^} we find 
2?A = 2? as in (|4.1ip . and that A is scaled stable if and only if 'jl/\{A) < 1. 

The state-space _ff ""-problem. The problems of finding tractable necessary and 
sufficient conditions for the strict state-space iJ°°-problem are similar to that for 
the state-space stabilizability problem. Here one also typically resorts to a more 
conservative 'scaled' version of the problem. 

We say that the realization {A, B, C, D} with decomposition (|4.8|) has scaled 
performance whenever there exists an invertible Q G I? so that 

9-1 

Iz®y 

or, equivalently, if there exists an X > in 2? so that 
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B ' 




' Q ■ 
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Iw®u 





< 1, 



(4.24) 
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< 0. 



(4.25) 



The equivalence of the two definitions goes as for the scaled stability case through 
the relation X — QQ* . Looking at the left upper entry in (|4.25p it follows that 
scaled performance of {A, B, C, D} implies scaled stability. Moreover, if (|4.24p 
holds for Q € V, then it is not hard to see that the transfer function G{z) in (|4.5p 
is also given by 

G{z) + C'{I - Z{z)A')-^Z{z)B' 
where the system matrix 

Q-i 



" A' 


B' ' 
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D 











■ A 


B ' 




' Q " 




C 


D 




Iw®u 



is equal to a strict contraction. It then follows from a standard fact on feedback 
connections (see e.g. Corollary 1.3 page 434 of 62\ for a very general formulation) 

that ||G'(z)|| < 1 for z G U'', i.e., G has strict performance. The scaled H°° -problem 
is then to find a controller K with realization {Ak, Bk, Ck, Dk} so that the closed 
loop system {Aci,Bci, Gci,Dci} has scaled performance. The above analysis shows 
that solving the scaled iJ°°-problem implies solving that state-space 7J°°-problem. 
The converse is again not true in general. Further elaboration of the same tech- 
niques as used in the proof of Theorem 14. 51 vields the following result for the scaled 
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7f°°-problem; see [THUSS]- For the connections between the Theorems 14.81 and 14.51 
in the more general setting of LFT models with structured uncertainty, we refer to 
|25j . Note that the result collapses to Theorem 12. 51 given in the Introduction when 
we specialize to the single- variable case d = 1. 

Theorem 4.8. Let {A, B, C, D} be a given realization. Then there exists a solution 
for the scaled H°° -problem associated with {A, _B,C, I?} if and only if there exist 
X,Y (ET> satisfying LMIs: 





No 




AY A* - Y 
CiYA* 
Bl 

'A*XA - X 
BIXA 



and the coupling condition 



AYCt 
CiYCl - 
n* 

A*XBi 
BlXBi - 



X 
I 



I 



Bl' 

Dii 
-I 

CI 
Dt 



I 
Y 



> 0. 





No 




< 0, r > 0, 



<o, x>o. 



(4.26) 



(4.27) 



(4.28) 



Here Nc and No are matrices chosen so that 

Nc is injective and Im Nc — Ker [Sj 
No is injective and Im No — Ker [C2 



D21] 



and 



Note that Theorem 14.81 docs not require that the problem be first brought into 
model-matching form; thus this solution bypasses the Nevanlinna-Pick-interpolation 
interpretation of the _ff°°-problem. 



4.3. Equivalence of frequency-domain and state-space formulations. In 

this subsection we suppose that we are given a transfer matrix G of size {nz + 
ny) X (nyv + nu) with coefficients in Q{C{z)ss) as in Section 14.11 with a given 
state-space realization as in Subsection 14.21 



G{z) = 



Gil 


G12 




'Dii 


D12 




'ci 


G21 


G22 




D21 


D22_ 


+ 





{I - Z{z)A)-^Z{z)[Bi B2] (4.29) 



where Z{z) is as in (|4.6p . We again consider the problem of finding stabilizing 
controllers K , also equipped with a state-space realization 



K{z) = Dk + Ck{I - Zk{z)AkT^Zk{z)Bk, 



(4.30) 



in either the state-space stability or in the frequency-domain stability sense. A 
natural question is whether the frequency-domain iJ°°-problem with formulation 
in state-space coordinates is the same as the state-space _ff°°-problem formulated 
in Section [321 

For simplicity in the computations to follow, we shall always assume that the 
plant G has been normalized so that D22 = 0. In one direction the result is clear. 
Suppose that K{z) = Dk + Ck{I — Z{z)Ak)~^ Z[z)Bk is a stabilizing controller 
for G{z) in the state-space sense. It follows that the closed-loop state matrix 



Aci = 



A 



- B2DKC2 

BkG2 
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B2CK 

Ak 



(4.31) 



is stable, i.e., / — Zci{z)Aci is invertible for all z in the closed polydisk D , with 
Zci{z) as defined in Subsection 14.21 On the other hand one can compute that the 



transfer matrix Q{G22,K) 



I -K(z 

-G22(z) I 



has realization 



W[z) 



'I Id,,' 




DkC2 Ck 
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C2 



Zcl{z)Acl) Zci(z) 



B2 B2DK 
Bk 



(4.32) 



As the resolvent expression {I — Zci{z)Aci) has no singularities in the closed 

polydisk D'', it is clear that W{z) has matrix entries in C(z)ss, and it follows that 
K stabilizes G22 in the frequency-domain sense. Under the assumption that G is 
internally stabilizable (frequency-domain sense) , it follows from Corollary 13.41 that 
K also stabilizes G (frequency-domain sense). 

We show that the converse direction holds under an additional assumption. The 
early paper [88] of Kung-Levy-Morf-Kailath introduced the notion of modal control- 
lability and modal observability for 2-D systems. We extend these notions to N-Y) 
systems as follows. Given a Givone-Roesser output pair {C, ^4}, we say that {C, A} 
is modally observable if the block-column matrix [ ^-■^^)^ j has maximal rank nx for 
a generic point z on each irreducible component of the variety det(/ — Z{z)A) — 0. 
Similarly we say that the Givone-Roesser input pair {A, i?} is modally controllable 
if the block-row matrix [i~AZ(z) b] has maximal rank nx for a generic point on 
each irreducible component of the variety det(/ — AZ{z)) = det(/ — Z{z)A) = 0. 
Then the authors of [88] define the realization {A, B, C, D} to be minimal if both 
{C, A} is modally observable and {A, B} is modally controllable. While this is a 
natural notion of minimality, unfortunately it is not clear that an arbitrary realiza- 
tion {A, B, C, D} of a given transfer function S{z) = D + C{I - Z{z)A)-^ Z{z)B can 
be reduced to a minimal realization {Aq, Bq^Cq, Dq} of the same transfer function 
Siz) = Do + Co{I - Z{z)Ao)-^Z{z)Bo. 

As a natural modification of the notions of modally observable and modally 
controllable, we now introduce the notions of modally detectable and modally sta- 
bilizable as follows. For {C, A} a Givone-Roesser output pair, we say that {C, ^} 
is modally detectable if the column matrix has maximal rank nx for a 

generic point z on each irreducible component of the variety det(/ — Z{z)A) — 

which enters into the polydisk D . Similarly, we say that the Givone-Roesser in- 
put pair {A, B} is modally stabilizable if the row matrix [i-AZ{z) b] has maxi- 
mal rank nx for a generic point z on each irreducible component of the variety 
det(/ — Z{z)A) = which has nonzero intersection with the closed polydisk D . 
We then have the following partial converse of the observation made above that 
state-space internal stabilization implies frequency-domain internal stabilization; 
this is an A^-D version of Theorem 12.61 in the Introduction. 



Theorem 4.9. Let and K3U\i be given realizations for G : [ ^ ] — > [ ^ ] 

and K: y ~*IA. Assume that {C2,A} and {Ck,Ak} are modally detectable and 
{A, B2} and {Ak, Bk} are modally stabilizable. Then K internally stabilizes G22 in 
the state-space sense ( and thus state-space stabilizes G) if and only if K stabilizes 
G22 in the frequency- domain sense (and G if G is stabilizable in the frequency- 
domain sense). 
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Remark 4.10. As it is not clear that a given realization can be reduced to a 
modally observable and modally controllable realization for a given transfer func- 
tion, it is equally not clear whether a given transfer function has a modally de- 
tectable and modally stabilizable realization. However, in the case that d — I, 
such realizations always exists and Theorem 14.91 recovers the standard 1-D result 
(Theorem 12.61 in the Introduction). 

The proof of Theorem 14.91 will make frequent use of the following basic result 
from the theory of holomorphic functions in several complex variables. For the 
proof we refer to [1281 Theorem 4 page 176]; note that if the number of variables d 
is 1, then the only analytic set of codimension at least 2 is the empty set and the 
theorem is vacuous; the theorem has content only when the number of variables is 
at least 2. 

Theorem 4.11. Principle of Removal of Singularities Suppose that the complex- 
valued function (p is holomorphic on a set S contained in of the form S — D — £ 
where T> is an open set in and £ is the intersection with T) of an analytic set of 
codimension at least 2. Then has analytic continuation to a function holomorphic 
on all of v. 

We shall also need some preliminary lemmas. 

Lemma 4.12. (1) Modal detectability is invariant under output injection, i.e., 
given a Givone-Roesser output pair {C, A} (where A . X ^ X and C : X ^ 
y) together with an output injection operator L: y X, then the pair 
{C, A} is modally detectable if and only if the pair {C, A + LC} is modally 
detectable. 

(2) Modal stabilizability is invariant under state feedback, i.e., given a Givone- 
Roesser input pair {A, _B} (where A: X ^ X and B.IA — > X) together 
with a state-feedback operator F : X ^ lA, then the pair {A, _B} is modally 
stabilizable if and only if the pair {A -\- BF, B} is modally stabilizable. 

Proof. To prove the first statement, note the identity 



7 -Z{z)L 




7 - Z{z)A 




I - 


Z{z){A + LC) 


/ 




C 






C 



Since the factor j ig invertible for all z, we conclude that, for each z G C^, 

[^~^^^^] maximal rank exactly when \^i-z{z)(A+lc)^^ Yyqs maximal rank, and 
hence, in particular, the modal detectability for {C, A} holds exactly when modal 
detectability for {C,A + LC) holds. 

The second statement follows in a similar way from the identity 

= [/ - (A + BF)Z{z) B] . 

□ 

Lemma 4.13. Suppose that the function W{z) is stable (i.e., all matrix entries of 
W are in C{z)ss) and suppose that 

W{z) = D + C(I - Z{z)A)-^Z{z)B (4.33) 

is a realization for W which is both modally detectable and modally stabilizable. 
Then the matrix A is stable, i.e., (/ — Z{z)A)~^ exists for all z in the closed 
polydisk ID) . 
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[/ - AZ{z) B] 



Proof. As W is stable and Z{z)B is trivially stable, then certainly 



■/ - Z{z)A 
C 



[I ^ Z{z)A)-^Z{z)B 



Z(z)B 
W{z) - D 



(4.34) 



is stable (i.e., holomorphic on J^f"). Trivially '^cf'^^\ maximal rank nx for 

all z G D*^ where det(/ — Z{z)A) ^ 0. By assumption, has maximal rank 

generically on each irreducible component of the zero variety of det(/ — Z{z)A) 
which intersects D . We conclude that ■^^■'^j has maximal rank nx at all 

points of D'' except those in an exceptional set £ which is contained in a subvariety, 
each irreducible component of which has codimension at least 2. In a neighborhood 
of each such point z e D — £, [ ^~-^^)'^] has a holomorphic left inverse; combining 
this fact with the identity (I4.34p . we see that (/ — Z{z)A)^^ Z{z)B is holomorphic 
on - 8. By Theorem iHl it follows that (/ - Z{z)A)-^ Z{z)B has analytic 

continuation to all of D . 
We next note the identity 

[Z(z) {I - Z{z)A)-^Z{z)B]= Z{z){I - AZ{z))-^[l - AZ{z) B] (4.35) 

where the quantity on the left-hand side is holomorphic on D by the result es- 
tablished above. By assumption {A, S} is modally stabilizable; by an argument 
analogous to that used above for the modally detectable pair {C, A}, we see that 
the pencil [/ — AZ{z) _B] has a holomorphic right inverse in the neighborhood of 

each point z in D — £' where the exception set £' is contained in a subvariety each 
irreducible component of which has codimension at least 2. Multiplication of the 
identity (|4.35[) on the right by this right inverse then tells us that Z{z){I — Z{z)A)~^ 

is holomorphic on D'' — £' . Again by Theorem 14.111 we conclude that in fact 
Z{z){I - Z{z)A)-'^ is holomorphic on all of D"*. 



We show that (/ — Z{z)A) is holomorphic on D as follows. Let Ej : X ^ Xj 
be the projection on the j-th component oi X = Xi ® ■ ■ ■ (B Xd- Note that the first 
block row of (/ — Z{z)A)~^ is equal to ziEi{I — Z{z)A)^^. This is holomorphic 

on the closed polydisk D*^. For z in a sufficiently small poly disk |zi| < p for i = 
l,...,(i, (/ — Z{z)A)~^ is analytic and hence ziEi{I — Z(z)A)~^ |zj=o = 0. By 
analytic continuation, it then must hold that zi{Ei{I — Z{z)A)~^ — for all z — 
(0, Z2, ■ • ■ , Zd) with \zi\ < 1 for i = 2, . . . ,d. For each fixed (z2, . . . , z^), we may 
use the single- variable result that one can divide out zeros to conclude that Ei {I — 
Z{z)A)~^ is holomorphic in zi at zi = 0. As the result is obvious for zi ^ 0, we 

conclude that Ei{I — Z{z)A)^^ is holomorphic on the whole closed polydisk Ui^. 
In a similar way working with the variable Zi, one can show that Ei{I — Z{z)A)~^ 
is holomorphic on the whole closed polydisk, and it follows that (/ — Z{z)A)^^ ~ 



(/ — Z{z)A) ^ is holomorphic on the whole closed polydisk as wanted. □ 



We are now ready for the proof of Theorem 

Proof of Theorem \4.9\ Suppose that K stabilizes G22 in the frequency-domain sense. 
This simply means that the transfer function W given by (|4.32p is holomorphic on 
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the closed polydisk D*^ . To show that Ad is stable, by Lemma 14.131 it suffices to 
show that {[-"gf' '^(f ] ,Aci} is modally detectable and that {Ad, is 
modally stabilizable. 

To prove that { [ ^^^^ ] ,Ad} is modally detectable, from the definition (|4?T7l) 
we note that 



A, 



'A ■ 




'B2 ■ 




DkC2 Ck 


Ak 


+ 


Bk 




C2 



By Lemma [4. 121 we see that modal detectability of { [ ' 



, Ad } is equivalent 



to modal detectability of { [ ] ' [ o t ] } • As [ 4' ] = [ a\ ] 

with [^j^ p] invertible, it is easily seen that modal detectability of the input pair 
{ [ ] ' [ A ] } is equivalent to modal detectability of { S^^Ait])- 

But the modal detectability of this last pair in turn follows from its diagonal form 

and the assumed modal detectability of {C2, A] and {Ck, Ak}- 

' B2 B2DK " 
.0 Bk . 



The modal stabilizability of {Ad, 
use of the identities 



I } follows in a similar way by making 



'A ■ 




B2DK B2 




'C2 


■ 




'B2DK B2 




'B2 ■ 




'Dk I 


Ak 


+ 


Bk 







Ck_ 


7 


Bk 




Bk 




I 



Ad = 
and noting that 



' Dk n 
. I oJ 



is invertible. 



□ 



In both the frequency-domain setting of Section 14.11 and the state-space set- 
ting of Section 14. 2( the true 7J°°-problem is intractable and we resorted to some 
compromise: the Schur-Agler-class reformulation in Section 14.11 and the scaled- 
/f°°-problem reformulation in Section 14.21 We would now like to compare these 
compromises for the setting where they both apply, namely, where we are given 
both the transfer function G and the state-space representation {A, B, C, D} for 
the plant. 

Theorem 4.14. Suppose that G(z) — [qH '^J^ ] is in model-matching form with 
state-space realization G{z) = D + C{I — Z{z)A)^^ Z{z)B as in (j4.29p . Suppose 
that the controller K{z) = Dk + Ck{I — Zk{z)Ak)^'^ Zk{z)Bk solves the scaled 
H°° -problem. Then the transfer function W{z) as in (j4.32p is a Schur-Agler-class 
solution of the Model-Matching problem. 

Proof. Simply note that, under the assumptions of the theorem, W{z) has a real- 
ization W = Dd -\- Cd{I — Zd{z)Ad)~^ Zd{z)Bd for which there is a state-space 
change of coordinates Q G P transforming the realization to a contraction: 



'A' 
C' 



B'' 
D 



< 1 where 



A' 
G' 



B'' 
D 







Ad 
Cd 



Bd 
Dd 







Thus we also have W{z) = D + G'{I - Zd{z)A')Zd{z)B' from which it follows 
that W is in the strict Schur-Agler class, i.e., < 1 for any d-tuple X — 

{Xi,...,Xd) of contraction operators Xj on a separable Hilbert space X. By 
construction W necessarily has the model matching form W = Gn -f G12AG21 
with A stable. □ 

Remark 4.15. In general a Schur-Agler function S{z) can be realized with a 
colligation matrix [ ^ £ ] which is not of the form 

(4.36) 



'A 


B 






0' 




'A' 


B'' 




Q 


0" 


C 


D 







/ 




G' 


D 







/ 
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with [ ^' ] equal to a strict contraction and Q £ P invertible. As an example, 
let A be the block 2x2 matrix given by Anderson-et-al in [1^. This matrix has 

2 

the property that / — Z{z)A is invertible for all z S D , but there is no Q S 2? so 

that ||Q~^A(5|| < 1. Here Z{z) and V are compatible with the block decomposition 

of A. Then for 7 > sufficiently small the function S{z) = 7(1 — Z{z)A)~^ has 

2 

||'S'(z)|| < p < 1 for some < p < 1 and all z £ D . Hence 5 is a strict Schur-class 
function. As mentioned in Section [4. li a consequence of the Ando dilation theorem 
[T7J is that the Schur class and the Schur-Agler class coincide for d = 2; it is not 
hard to see that this equality carries over to the strict versions and hence S is in 
the strict Schur-Agler class. As a consequence of the strict Bounded-Real-Lemma 
in [53], S admits a strictly contractive state-space realization [p/ ^ ]• However, 
the realization [^^ Bl ~ [7/7/] obtained from the fact that 

S{z) - 7(/ - Z{z)A)-' = 7/ + 7(/ - Z{z)A)-'Z{z)A, 

cannot relate to [ ^ ] as in (|4.36p since that would imply the existence of an 
invertible Q G I? so that Q'^AQ = A' is a. strict contraction. 

Remark 4.16. Let us assume that the G{z) in Thcorem l4.14l is such that G12 and 
G21 are square and invertible on the distinguished boundary T"^ of the polydisk D"* 
so that the Model-Matching problem can be converted to a polydisk bitangential 
Nevanlinna-Pick interpolation problem along a subvariety as in [32j . As we have 
seen, the solution criterion using the Agicr interpolation theorem of [11 [35] then 
involves an LOI (Linear Operator Inequality or infinite LMI). On the other hand, 
if we assume that we are given a stable state-space realization {A, B, C, D} for 

G{z) = Cailz) ^^0^^^ , we may instead solve the associated scaled iJ°°-problem 
associated with this realization data-set. The associated solution criterion in The- 
orem |TS] remarkably involves only finite LMIs. A disadvantage of this state-space 
approach, however, is that in principle one would have to sweep all possible (sim- 
ilarity equivalence classes of) realizations of G{z); while each non-equivalent re- 
alization gives a distinct i/°°-problem, the associated frequency-domain Model- 
Matching/bitangential variety-interpolation problem remains the same. 

4.4. Notes. In [S^ Lin conjectured the result stated in Theorem 14.11 that G22- 
stabilizability is equivalent to the existence of a stable coprime factorization for 
6*22- This conjecture was settled by Quadrat (see [122[ 11171 1120] ) who obtained 
the equivalence of this property with projective-freeness of the underlying ring and 
noticed the applicability of the results from [46l [83] concerning the projective- 
freeness of C(z)ss- 

For the general theory of the A^-D systems, in particular for N=2, considered in 
Subsection [421 we refer to [811155] . 

The sufficiency of scaled stability for asymptotic/Hautus-stability goes back to 
[5S] . Theorem l4.5l was proved in ^5] for the more general LFT models in the context 
of robust control with structured uncertainty. The proof given here is based on the 
extended Finsler's lemma fLemma I4.7p . and basically follows the proof from [S5] 
for the solution to the scaled i/°°-problem (Theorem 14. 8p . As pointed out in [5S] . 
one of the advantages of the LMI-approach to the state-space problem, even 
in the classical setting, is that it allows one to seek controllers that solve the scaled 
iJ°°-problem with a given maximal order. Indeed, it is shown in [66l [18] (see also 
[57] ) that certain additional rank constraints on the solutions X and Y of the LMIs 
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(|4.26p and (|4.27p enforce the existence of a solution with a prescribed maximal 
order. However, these additional constraints destroy the convexity of the solution 
criteria, and are therefore usually not considered as a desirable addition. 

An important point in the application of Finsler's lemma in the derivation of the 
LMI solution criteria in Theorems 14.51 and 14 . 81 is that the closed- loop system matrix 
Aci in (|4.3ip has an affine expression in terms of the unknown design parameters 
{Ak, Bk, Ck, Dk}- This is the key point where the assumption D22 = is used. A 
parallel simplification occurs in the frequency-domain setting where the assumption 
G22 — leads to the Model-Matching form. The distinction however is that the 
assumption G22 = is considered unattractive from a physical point of view while 
the parallel state-space assumption D22 G22(0) = is considered innocuous. 

There is a whole array of lemmas of Finsler type; we have only mentioned the 
form most suitable for our application. It turns out that these various Finsler lem- 
mas are closely connected with the theory of plus operators and Pesonen operators 
on an indefinite inner product space (see [44]). An engaging historical survey on 
all the Finsler's lemmas is the paper of Uhlig [135j . 

The notions of modally detectable and modally stabilizable introduced in Sub- 
section 14.31 along with Theorem 14.91 seem new, though of somewhat limited use 
because it is not known if every realization can be reduced to a modally detectable 
and modally stabilizable realization. We included the result as an illustration of 
the difficulties with realization theory for iV-D transfer functions. 

We note that the usual proof of Lemma 14.131 for the classical 1-D case uses 
the pole-shifting characterization of stabilizability/detectability (see [S71 Exercise 
2.19]). The proof here using the Hautus characterization of stabilizability/detectability 
provides a different proof for the 1-D case. 



5. Robust control with structured uncertainty: the commutative 

CASE 

In the analysis of 1-D control systems, an issue is the uncertainty in the plant 
parameters. As a control goal, one wants the control to achieve internal stability 
(and perhaps also performance) not only for the nominal plant G but also for a 
whole prescribed family of plants containing the nominal plant G. 

A question then is whether the controller can or cannot have (online) access to 
the uncertainty parameters. In a state-space context it is possible to find sufficient 
conditions for the case that the controller cannot access the uncertainty parameters, 
with criteria that are similar to those found in Theorems 14.51 and 14.81 but additional 
rank constraints need to be imposed as well, which destroys the convex character 
of the solution criterion. The case where the controller can have access to the 
uncertainty parameters is usually given the interpretation of gain-scheduling, and 
fits better with the multidimensional system problems discussed in Section ID In 
this section we discuss three formulations of 1-D control systems with uncertainty 
in the plant parameters, two of which can be given gain-scheduling interpretation, 
i.e., the controller has access to the uncertainty parameters, and one where the 
controller is not allowed to use the uncertainty parameters. 
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5.1. Gain-scheduling in state-space coordinates. Following [106j . we suppose 
that we are given a standard linear time- invariant input/state/output system 




= AM{Su)xit) 
= CMi{Su)xit)- 

= CM2{Su)xit)- 



but where the system matrix 



Cm2{Su) 



Bmi{5u) 
Dmii{&u) 
Dm2i{5u) 



rBMl{Su)w{t) ' 

DMiiiSu)w(t) 

DM2liSu)w(t) 



Bm2{Su) 

Dmi2{5u) 

Dm22{Su) 



BM2{5u)u[t) 

- DM12{Su)u{t) 

- DM22{Su)u{t) 



{t e Z+) 

(5.1) 













w 




z 




u 




y 



G{6) 



is not known exactly but depends on some uncertainty parameters 8u 
in C^. Here the quantities 8i are viewed as uncertain parameters which the con- 
troller can measure and use in real time. The goal is to design a controller Sa' 
(independent of (5c/) off-line so that the closed- loop system (with the controller ac- 
cessing the current values of the varying parameters <5i , . . . , (5^ as well as the value of 
the measurement signal y from the plant) has desirable properties for all admissible 
values of bu , usually normalized to be < 1 for k — 1, . . . , d. 

The transfer function for the uncertainty parameter 5u can be expressed as 

-Da/11 ((5[/) Dmi2{Su) 
_Dm2i{Su) Dm22{Su)_ 

Cmi{5u) 
_Cm2{Su)_ 

where we have introduced the aggregate variable 

5 = ((5(7, A) = ((5i, . . . , (5d, A). 

It is not too much of a restriction to assume in addition that the functional 
dependence on 5u is given by a linear fractional map (where the subscript U suggests 
uncertainty and the subscript S suggests shift) 

Ass 
Cis 
C2S 



+ \ 



{Ix ~ \AM{5u)r^ [Bmi{5u) Bm2{5u)\ (5.2) 



Am{5u) 
Cmi{Su) 

Cm2{Su) 



Bmi{Su) 
Dmii{5u) 
Dm2i{Su) 



Bm2{5u) 
Dmi2{5u) 

Dm22{5u) 



Bsi 
D21 



BS2 

D12 
D22 



Asu 
C2U 



{I - Z{5u)Auu)'^Z{5u)[ Aus Bui Bu2 ] 



where Z{Su) is defined analogously to Z{z) in (|4.6[) relative to a given decomposition 
of the "uncertainty" state-space Xjj — Xjj^i®- ■ -(BXij^d on which that state operator 
Auu acts. In that case the transfer function G{6) admits a state-space realization 

G{S)=B' ?'l = [n" + ^2] (5.3) 

[Lt21 (-T22J [i^21 ^22} L 2j 

with system matrix given by 



A 

Ci 
C2 



Bi 
D21 



B2 
D12 
D22 



' Auu 


Aus 


Bui 


Bu2 


Asu 


Ass 


Bsi 


BS2 


Ciu 


Cis 


Dn 


D12 


C2U 


C2S 


D21 


D22 



(5.4) 
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Here Z{6) is again defined analogously to (|4.6p but now on the extended state-space 
Xext = Xij ® X. 

We can then consider this gain-scheduling problem as a problem of the con- 
structed iV-D system (with N = d + 1), and seek for a controller K with a state- 
space realization 

K{d) = Dk + Ck{I - Zk{5)Ak)Zk{5)Bk (5.5) 

so that the closed loop system has desirable properties from a gain-scheduling per- 
spective. Making a similar decomposition of the system matrix for the controller 
K as in (|5.4p . we note that K{5) can also be written as 

where ^j\/,_f£-(5c/), Bm,k{5u), CM.Ki^u) and Dm.k^^u) appear as the transfer func- 
tions of A^-D systems (with N = d), that is, K{S) can be seen as the transfer 
function of a linear time-invariant input/state/output system 



XK{t+l) = AM,K{Su)xK{t) + BM.K{5u)u{t) , „ 

u{t) = CM,K{5u)xK{t) + DMM5u)y{t) ^"^^ 



depending on the same uncertainty parameters 5jj = (^i, . . . , 5d) as the system E. 

Similarly, the transfer function Gci{5) of the closed- loop system with system 
matrix [ '^"'■^ ^"'■^ ] as defined in (|4.10p also can be written as a transfer matrix 

Gci{5) - DmA^u) + \CmA5u){I - XAMA5u)r^BMA5u) 

with Am,ci{^u), BmA^u)-, Cm,ci{Su) and Dm,ci{^u) transfer functions of A^-D sys- 
tems (with N — d), and the corresponding linear time- invariant input/state/output 
system 



x{t + l) = AMA^u)xit) + BM,cl{Mw{t) 
z{t) = CMA^u)x{t)+DMA^u)w{t) 



{n G Z+) 



also appears as the closed- loop system of E and Ex- 
it then turns out that stability of Ac/, that is, / — Zci{5)Aci invertible for all S in 

D ^ (with Zci as defined in Subsection l4.2|l corresponds precisely to robust stability 
of Ec/, i.e., the spectral radius of AM,ciiSu) is less than 1 for all Sjj — (6i, . . . , Sd) 
so that \Sk\ < 1 for k = 1, . . . ,d, and K with realization (|5.5p solves the state-space 
i?°°-problem for G with realization (|5.3p means that the closed loop system Ec/ 
has robust performance, i.e., Ec/ is robustly stable and the transfer function Gd 
satisfies 

||Gc/(5)|| < 1 for aU 5 = (^1, . . . , 5^, A) e D''^\ 

We may thus see the state-space formulation of the gain-scheduling problems con- 
sidered in this subsection as a special case of the N-D system stabilization and 
i/°°-problems of Subsection 14.21 In particular, the sufficiency analysis given there, 
and the results of Theorem 14.51 and 14.81 provide practical methods for obtaining 
solutions. As the conditions are only sufficient, solutions obtained in principle may 
be conservative. 
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5.2. Gain-scheduling: a pure frequency-domain formulation. In the ap- 
proach of Helton (see [731 [H]), one eschews transfer functions and state-space co- 
ordinates completely and supposes that one is given a plant G whose frequency 
response depends on a load with frequency function 5{z) at the discretion of the 
user; when the load 5 is loaded onto G, the resulting frequency-response function 
has the form G{z,5{z)) where G = G(-, •) is a function of two variables. The control 
problem (for the company selling this device G to a user) is to design the controller 
K = K{-,-) so that K{-,5{-)) solves the iJ°°-problem for the plant G(-,(5(-)). The 
idea here is that once the user loads 5 onto G with known frequency-response func- 
tion, he is also to load 5 onto the controller K (designed off-line); in this way the 
same controller works for many customers using many different (5's. When the dust 
settles, this problem reduces to the frequency-domain problem posed in Section [4. II 
with d = 2; an application of the Youla-Kucera parametrization (or simply using 
the function Q{z) — K{z){I G 22(2) K{z))~^ if the plant G itself is stable) reduces 
the problem of designing the control K to a Nevanlinna-Pick-type interpolation 
problem on the bidisk. 



5.3. Robust control with a hybrid frequency-domain/state-space formu- 
lation. We now consider a hybrid frequcncy-domain/state-space formulation of the 
problem considered in Subsection 15.11 the main difference is that in this case the 
controller is not granted access to the uncertainty parameters. 

Assume we are given a 1-D-plant G(A) that depends on uncertainty parameters 
Su = (Si, ■ ■ ■ , Sd) via the linear fractional representation 



G{Su,X) 



Gii(A) Gi2(A) 

G2i(A) G22(A) 

'GwW 
G2u{X) 



iI^ZiSu)GuuW)-'ZiSu)[Guii\) Gu2i\)] (5.6) 



with Z{Su) as defined in Subsection 15. H and where the coefficients are 1-D-plants 
independent of Su- 



Gaug (A) 



GuuW Gt/i(A) G[/2(A) 
Gn/(A) Gii(A) Gi2(A) 
G2c/(A) G2i(A) G22(A) 







Xjj 


w 




Z 


u 




y 



In case GaugW is also given by a state-space realization, we can write G{du, A) as in 
(|5.3p with S = {Sjj, A) and Z{d) acting on the extended state-space X^xt = '^u © X. 

For this variation of the gain-scheduling problem we seek to design a controller 
K{X) with matrix values representing operators from 3^ to W so that K solves 
the _ff°°-problem for G{Su,X) for every 6u with ||Z((5c/)|| < 1, i.e., \6j\ < 1 for 
j ~ 1, . . . , d. For the sequel it is convenient to assume that Z — W. In that case, 
using the Main Loop Theorem [141[ Theorem 11.7 page 284], it is easy to see that 
this problem can be reformulated as: Find a single-variable transfer matrix K(-) 



so that Q{G,K) given by ((2?2)) . with G = 



Gii(A) Gi2(A) 
G2i(A) G22(A) 



Gi 
G2 



2.21) taken to be 



Guu{X) 
Giu{X) 


Gui{\) 
Gii(A) 


Gu2{X) ' 
Gi2(A) 


. G2[/(A) 


G2i(A) 


G22(A) 
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is stable and such that 

AiA (Gii(A) + Gi2(A)(/ - K{X)G22{X)r'K{X)G2i{X)) < 1. 
Here //a is as defined in (|4.22p witli A the C* -algebra 

Application of the Youla-Kucera parametrization of the controllers K that sta- 
bihze 0(G,if) as in Subsection 13.31 converts the problem to the following: Given 
stable 1-variable transfer functions 7i(A), T2{X), and T^^X) with matrix values rep- 
resenting operators in the respective spaces 

find a stable 1-variable transfer function A(A) with matrix values representing op- 
erators in C{Xjj © Xjj © U) so that the transfer function S{X) given by 

5(A) = Ti(A) + r2(A)A(A)T3(A) (5.7) 

has n^{S{X)) < 1 for all A G D. If T2(C) and T3(C) are square and invertible for 
( on the boundary T of the unit disk D, the model-matching form (|5.7p can be 
converted to bitangential interpolation conditions (see e.g. [261 ) : for simplicity, say 
that these interpolation conditions have the form 

XtS{Xi)=yi, S{X'j)uj = Vj for i ^ I, ... ,k, j^l,...,k' (5.8) 

for given distinct points ^i,X'j in B, row vectors Xi,yi and column vectors Uj,Vj. 
Then the robust iJ°°-problcm {H°^ rather than rational version) can be converted 
to the /i-Nevanlinna-Pick problem: find holomorphic function S on the unit disk 
with matrix values representing operators in C{Xu W, Xjj © Z) satisfying the 
interpolation conditions (j5.8p such that also 

MA (5(A)) < 1 for all A G D. 

It is this /i-version of the Nevanlinna-Pick interpolation problem which has been 
studied from various points of view (including novel variants of the Commutant 
Lifting Theorem) by Bercovici-Foias-Tannenbaum (see [551 [5^ HUl l¥T| ) and Agler- 
Young (see [H [3 [HI E] and Huang-Marcantognini- Young [77j). These authors 
actually study only very special cases of the general control problem as formulated 
here; hence the results at this stage are not particularly practical for actual control 
applications. However this work has led to interesting new mathematics in a number 
of directions: we mention in particular the work of Agler- Young on new types of 
dilation theory and operator- model theory (see [51 [S]), new kinds of realization 
theorems [TU], the complex geometry of new kinds of domains in (see [51 [TH 
I13p. and a multivariable extension of the Bercovici-Foias-Tannenbaum spectral 
commutant lifting theorem due to Popescu [114| . 

5.4. Notes. In the usual formulation of /i (see [10711141] ). in addition to the scalar 
blocks 5iln in Z{5), it is standard to also allow some of the blocks to be full 

5« - 5<!?. 



blocks of the form A, 



The resulting transfer functions then 



1:1 ■■■ 

have domains equal to be (reducible) Cartan domains which are more general than 
the unit polydisk. The theory of the Schur-Agler class has been extended to this 
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setting in p31 120j . More generally, it is natural also to allow non-square blocks. A 
formalism for handling this is given in [29j ; for this setting one must work with the 
intertwining space of A rather than the commutant of A in the definition of in 
(|4.23p . With a formalism for such a non-square uncertainty structure available, one 
can avoid the awkward assumption in Subsection 15.31 and elsewhere that W ~ Z. 



6. Robust control with dynamic time- varying structured uncertainty 



6.1. The state-space LFT-model formulation. Following [271 [Ml [Ml [THH] , we 

now introduce a variation on the gain-scheduling problem discussed in Section 15.11 
where the uncertainty parameters djj — . . . ,6^) become operators on the 
space of square-summable sequences of complex numbers indexed by the integers 
Z, and are to be interpreted as dynamic, time-varying uncertainties. To make 
the ideas precise, we suppose that we are given a system matrix as in (I5.4|l . We 
then tensor all operators with the identity operator 1^2 on i'^ to obtain an enlarged 
system matrix 





'A Bi 


B2' 


M = 


Ci Dn 


D12 




C2 D21 


D22_ 


which we also write as 






A Bi 


B2 


M = 


Ci Dn 


D12 




C2 D21 


D22 





' Auu 


Aus 


Bui 


Bu2 




Asu 


Ass 


Bsi 


Bs2 




Ciu 


Cis 


Dn 


D12 




C2U 


C2S 


D21 


D22 



® Ip, 



'{Xu © Xs) ® l^' 



\Xu (B Xs) e^' 



(6.1) 



(6.2) 



Given a decomposition Xjj — Xjji © • • • © Xud of the uncertainty state space Xjj^ we 
define the matrix pencil Xij{5u) with argument equal to a d-tuple djj = (5i, . . . , 5d) 
of (not necessarily commuting) operators on P by 



Zc/(^t 



IXui ® ^1 



In addition we let S denote the bilateral shift operator on we sometimes will 
also view S as an operator on the space i of all sequences of complex numbers or 
on the subspace of l'^ that consists of all sequences in £^ with finite support. 
We obtain an uncertain linear system of the form 




where the system matrix 



Am{Su)x + Bmi{Su)w + Bm2{Su)u 
Cmi{Su)x 4- Dmii{Su)w + DMi2iSu)u 
Cm2{^u)x + Dm2i{Su)w + Dm22{Su)u 



(6.3) 



Am{Su) 


Bmi{^u) 


Bm2{^u) 










' Xs®i' 




Dmii{Su) 


Dm12{^u) 












Cm2{Su) 


Dm2i{Su) 


Dm22{^u)_ 






5^ 
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is obtained from the feedback connection 







xu 




xs 


= M 


XS 


, subject to Xu 


z 




z 


. y . 




y 





Zu{Su)xu, 



that is, 



Am(^c/) Bmi{Su) 



Cmi{Su) 


Dmii 


iSu) 


Cm2{^u) 


Da/21 










+ 


Clu 


(1- 




C2U 





Ass 


Bsi 


Bs2 




Cis 


Dn 


D12 


+ 


C2S 


D21 


D22 




(Su) 


[ A[/s 




Bc/2 



'Bm2{Su) 
^M12{Su) 
^M22{^u) 



ZuiSu)Auu)-'Zi 



(6.4) 

As this system is time-varying, due to the presence of the time-varying uncer- 
tainty parameters Su, it is not convenient to work with a transfer- function acting 
on the frequency-domain; instead we stay in the time-domain and work with the 
input-output operator which has the form 



Da/11 (^y) Da/12 (^t/) _^ Cmi{Su 

Gm2i{Su) Da/22(^C/)J [Ca/2(^;7 

X {Ixs»e^ - {Ixs ® S)Aa/(^[/))"^ 



(6.5) 

Ba/2(<5[/)] , 



Now write 8 for the collection {5u: S) of d -I- 1 operators on l"^ . Then the input- 
output operator G{8) given by (j6.5p has the noncommutative transfer-function 
realization 



G{8) = 





G12W 




Dn 


D12 




Ci 


G2l(^) 


G22(^) 




D21 


D22 


+ 





Z[/(<Sc/) 
/. 



(/-Z(,5)A)-iZ(J) [Bi B2] 

(6.6) 

In the formulas (16.411- 



with system matrix as in (|6.ip and Z(5) 

()6.6p the inverses may have to be interpreted as the algebraic inverses of the corre- 
sponding infinite block matrices; in that way, the formulas make sense at least for 
the nominal plant, i.e., with 5u = (0, . . . , 0). 

More generally, the transfer-function G can be extended to a function of d + 1 
variables in C{1'^) by replacing S with another variable 5d+i G 'C(^^). In that case, 
the transfer-function can be viewed as an LFT-model with structured uncertainty, 
as studied in [OlllS^. However, as a consequence of the Sz.-Nagy dilation theory, 
without loss of generality it is possible in this setting of LFT-models to fix one of the 
variables to be the shift operator S; in this way the LFT-model results developed 
for d+1 free variable contractions apply equally well to the case of interest where 
one of the variables is fixed to be the shift operator. 

Such an input/state/output system S with structured dynamic time- varying 
uncertainty 5u is said to be robustly stable (with respect to the dynamic time- 
varying uncertainty structure TiuiSu)) if the state-matrix Aa/(^[/) is stable for all 
choices of 5u subject to ||Zc/(^//)|| < 1, that is, if Ixsm'^ ~ i^Xs ^ S)Aa/(^(7) is 
invertible as an operator on Xs © for all du with ||Z(7(5(7)|| < 1. Since 



AuiSu) = Ass + Asc/(/ - ZuiSu)Auur'Zidu)A 



USi 
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it follows from the Main Loop Theorem [141i, Theorem 11.7 page 284], that this 
condition in turn reduces to: 



Ix - Z{S)A is invertible for all S = (<5[/, S) with \\Z{6)\\ < 1. 



(6.7) 



Note that this condition amounts to a noncommutative version of the Hautus- 
stability criterion for the matrix A (where A = A ) . We shall therefore call the 
state matrix A nc-Hautus-stable if (j6.7p is satisfied (with nc indicating that we are 
in the noncommutative setting). The input/state/output system S is said to have 
nc-performance (with respect to the dynamic time-varying uncertainty structure 
Zui^u)) if it is robustly stable (with respect to this dynamic time-varying uncer- 
tainty structure) and in addition the input-output operator G{d) has norm strictly 
less than 1 for all choices of ^ = {djj, S) with ||Z(^)|| < 1. 

One of the key results from the thesis of Paganini [108J which makes the non- 
commutative setting of this section more in line with the 1-D case is that, contrary 
to what is the case in Subsection 14.21 for operators A = A Q) 1^2 on X ® i'^ we do 
have /iA(A) = 'p,^{A) when we take A to be the C*-algebra 











Xs 



>d+l 



(^1 



1 



(6.8) 

Write X> for the commutant of A in C{{Xu ® Xs)®i^)- Then the main implication 
of the fact that ii/^{A) ~ 'jl^{A) is that nc-Hautus-stability of A is now the 
same as the existence of an invertible operator Q G I? so that ||Q^^AQ|| < 1 or, 
equivalently, the existence of a solution X G I? to the LMIs A*XA — A < and 
X > 0. However, it is not hard to see that X is an element of X> if and only if 
X = X ® with X being an element of the C*-algebra V in (|4.1ip . Thus, in fact, 
we find that A — A® If2 is nc-Hautus-stable precisely when A is scaled stable, i.e., 
when there exists a solution X eV to the LMIs A*XA - A <0 and X > 0. 

These observations can also be seen as a special case (when C2 — and B2 = 0) 
of the following complete analogue of Theorem [231 for this noncommutative setting 
due to Paganini |108j . 

Proposition 6.1. Given a system matrix as in (|6.ip - (|6.2p . then: 

(i) The output pair {C2, A} is nc-Hautus-detectable, that is, for every 8 = 
{di,...,Sd+i), with 5j e C{£^) forj^l,...,d+l,so that \\Z{S)\\ < 1 the 
operator 



I 



Z(<5)A 
C2 



{Xu ® Xs) ® f 



{Xu ® Xs) » 



has a left inverse, if and only if {C2, A} is nc- operator- detectable, i.e., 
there exists an operator Ij = L ® Ip , with L : y —>■ X, so that A -\- LC2 is 
nc-Hautus-stable, if and only if there exists a solution X ^ V to the LMIs 



A*XA ~X - CW2 < 0, X>{). 



(6.9) 



(ii) The input pair {A,B2} is nc-Hautus-stabilizable, that is, for every S = 
{di,...,Sd+i), with e C{f) forj^l,...,d+l, so that \\Z{S)\\ < 1 the 
operator 

{Xu ® Xs) (E> e 
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[ I-Z{5)A B2 ] : 



[Xu ® Xs) «) 



has a left inverse, if and only i/{A, B2} is nc- operator- stabilizable, i.e., 
there exists an operator Y — F ® Ip , with F : X -^lA, so that A + B2F is 
nc-Hautus-stable, which happens if and only if there exists a solution Y £ V 
to the LMIs 



AY A* - Y - B2B; < 0, y > 0. 



(6.10) 



In case the input/state/output system S is not stable and/or does not have 
performance, we want to remedy this by means of a feedback with a controller 
K, which we assume has on-line access to the structured dynamic time- varying 
uncertainty operators du in addition to being dynamic, i.e., K ~ K(5) = K((5[/, S). 
More specifically, we shall restrict to controllers of the form 



K{d) = Da + Ck{I - ZKiS)AK)-'ZKiS)BK 



(6.11) 



where 







'-'Kui^U 



Ixki ® ^1 



IXk 



with system matrix Mk of the form 



'Ak 













(6.12) 



where Xku = ^Kui 
tensor-factorization 



{Xku © Xks) «) I 
S XkucLi and where the matrix entries in turn have a 





Ba 




'Ak 


® Ip Bk ® Ip 


Ca 


Da 




Ck 


Dk ® 



(6.13) 



If such a controller K(<5) is put in feedback connection with G(^), where we 
impose the usual assumption D22 = to guarantee well-posedness, the resulting 
closed-loop system input-output operator Gci{S), as a function of the operator 
uncertainty parameters 6u = (<5i, . . . , 8d) 
is formally exactly as in (|4.9p . that is 

with system matrix 



and the shift S, has a realization which 



Aci 



Be/ 



A 


■f B2DAC2 


B2CA' 


Bi 


+ B2DAD21 




BaC2 


Aa 




BaD21 


Ci 


f D12DAC2 


D12CA 


Dn 


+ D12DAD21 



(6.14) 



and 



Zc/(<5) 



where 5 = {6ij, S) 



(6.15) 



which is the same as the system matrix (I4.10p tensored with Ip , 

Z{8) 
Za(<5) 

The state-space nc- stabilization problem (with respect to the given dynamic time- 
varying uncertainty structure 5u) then is to design a controller K with state-space 
realization {Aa, Ba, Ck, Da} as above so that the closed-loop system Sci defined 
by the system matrix (I6.14p is robustly stable. The state-space nc-H°° -problem is 
to design a controller K with state-space realization {Aa:, Ba, Ca, Da} as above 
so that the closed-loop system Hci also has robust performance. 

Since the closed-loop state-operator Ad is equal to Ad ® Ip with Ad defined 
by (|4.10p . it follows as another implication of the fact that /^a is equal to /Ua for 

49 



operators that are tensored with Ip (with respect to the appropriate C*-algebra 
A) that Ac/ is nc-Hautus-stable precisely when Ad is scaled stable, i.e., we have 
the following result. 

Proposition 6.2. Let S andT, be the systems given by (j6.3p and (j5.ip . respectively, 
corresponding to a given system matrix (|5.4p . Then S is nc-Hautus-stabilizable if 
and only if S is scaled-stabilizable. 

Thus, remarkably, the solution criterion given in Section [4. 21 for the scaled state- 
space stabilization problem turns out to be necessary and sufficient for the solution 
of the dynamic time- varying structured-uncertainty version of the problem. 

Theorem 6.3. Let S be the system given by (16. 3p corresponding to a given sys- 
tem matrix (j6.f p . Then S is nc-Hautus-stabilizable if and only if the output pair 
{C2, A} is nc-Hautus- detectable and the input pair {A,B2} is nc-Hautus-stabiliz- 
able, i.e., if there exist solutions X,Y ^T>, with T> the C* -algebra given in (|4.f ip . 
to the LMLs and (|6l0l) . In this case K - [ ^^i^^ ] (g) Ip with [cId'k] ""^ 

in (|4.12p is a controller solving the nc-Hautus stabilization problem for S. 

In a similar way, the state-space nc-i/°°-problem corresponds to the scaled 77°°- 
problem of Subsection (|4.2p . 

Theorem 6.4. Let S be the system given by (j6.3p for a given system matrix (|6.ip . 
Then there exists a solution K, with realization ()6.1ip . to the state-space nc-H°°- 
problem for the non- commutative system S if and only if there exist X,Y ^T> that 
satisfy the LMIs (j4.27p and (j4.26p and the coupling condition (|4.28p . 

Proof. Let S and E be the systems given by (j6.3p and (|5.ip . respectively, corre- 
sponding to a given system matrix (|5.4p . Using the strict bounded real lemma from 
[29] in combination with similar arguments as used above for the nc-stabilizability 
problem, it follows that a transfer- function K with realization (|6.1ip - (|6.13p is a so- 
lution to the state-space nc-iJ°°-problem for S if and only if the transfer function 
K with realization (|4.7p is a solution to the scaled i?°°-problem for the system E. 
The statement then follows from Theorem 14.81 □ 

6.2. A noncommutative frequency-domain formulation. In this subsection 
we present a frequency-domain version of the noncommutative state-space setup of 
the previous subsection used to model linear input /state/output systems with LFT- 
model for dynamic time-varying structured uncertainty. The frequency-domain 

setup here is analogous to that of Section 14.11 but the unit polydisk D is replaced 
by the noncommutative polydisk consisting of all d-tuples 5 = {5i, . . . , 8d) of 
contraction operators on a fixed separable infinite-dimensional Hilbert space tC. 

We need a few preliminary definitions. We define !Fd to be the free semigroup 
consisting of all words a = in ■ ■ ■ ii in the letters {1, . . . , d}. When a = in ■ ■ ■ ii we 
write N — \a\ for the number of letters in the word a. The multiplication of two 
words is given by concatenation: 

a- (3 = iN--- iijM ■■■.]! if a = • ■ • ii and (5 = jm--- ji- 

The unit element of Td is the empty word denoted by with |0| = 0. In addition, 
we let z = (zi, . . . , Zd) stands for a d-tuple of noncommuting indeterminates, and 
for any a ~ iN ■ ■ ■ ii G .^d ~ {0}, we let z" denote the noncommutative monomial 
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= ■ ■ ■ Zi-^, while z^ — 1. If a and /3 are two words in J-d, we multiply the 
associated monomials and z^ in the natural way: 

z"-z^= z'^-P. 

Given two Hilbert spaces U and y, we let C{U,y){{z)) denote the collection of 
all noncommutative formal power series S{z) of the form S{z) = X^aej'^d SaZ°' 
where the coefficients Sa are operators in C{U,y) for each a £ Td- Given a 
formal power series S{z) = J2aeJ='a ^<^^" together with a d-tuple of linear operators 
S = . . . , 6d) acting on we define S{S) by 

S{5) = lim y Sa^^S'^ e C{U(E)IC,y(E> IC) 

aeJ^a - \a\=N 

whenever the limit exists in the operator-norm topology; here we use the notation 
5°' for the operator 

<5" = • • • J,, if a = ijv • ■ -ii e - {0} and - I,c- 

We define the noncommutative Schur-Agler class SAncd^l^jy) {strict noncom- 
mutative Schur-Agler class SA^cdi^^y)) to consist of all formal power series in 
C{U,y){{z)) such that \\S{d))\\ < 1 {\\S{6)\\ < 1) whenever <5 = {6i,...,Sd) is a 
d-tuple of operators on K- with \\5j || < 1 || < 1) for j = 1, . . . , d. Let 

B„,^d:={S = {S,,...,6d): 6^ e C{IC), \\S,\\<1, j = l,...,d}, 
Bnc,d-^{S^{Su...,Sd): 6jeC{IC), \\S,\\ < 1, j = 1, . . . , d}. 

We then define the strict noncommutative H°° -space H°°'°{C{U,y)) to consist of 
all functions F from D„c,rf to C{U ® K.,y ® K,) which can be expressed in the form 

m - s{5) 

for all 5 e Dnc.d where p^^S is in the strict noncommutative Schur-Agler class 
SA°„^^di^,y) for some real number p > 0. We write H^^^^{C{U,y)) for the set of 
functions G from Binc,d to C{U ® JC.y ® IC) that are also of the form G{5) = S{6), 
but now for S G ^nc,d and p~^S in SAnc,d{l^,y) for some p > 0. Note that 
SAnc,d{l^, y) amounts to SAnc,d{'C, C) (g) y). In the sequel we abbreviate the 
notation SAnc.di'C^C) for the scalar Schur-Agler class to simply SAnc,d- Similarly, 
we simply write 5^°,,^, and H^^ ^ instead of 5^°,.d(C, C), ff^'°(C,C) and 

^^.d(C,C), respectively. Thus we also have H^^l°{^{u]y)) = <^ -C(W,3^), 

etc. We shall be primarily interested in the strict versions SA^c d -^^ d '^^ the 
noncommutative Schur-Agler class and _ff°°-space. 

We have the following characterization of the space H^^'°^{C{U,y)). For the 
definition of completely positive kernel and more complete details, we refer to j30) . 
The formulation given here does not have the same form as in Theorem 3.6(2) of 
[30j , but one can use the techniques given there to convert to the form given in the 
following theorem. 

Theorem 6.5. The function F: JS>nc,d J0.{U (d IC^y ® K.) is in the strict non- 
commutative H°°-space H^^''^(C{U , y)) if and only if there are d strictly completely 
positive kernels 

Kk : (D„,,d X D„,,d) X C{JC) C{Y ^ IC) for k ^ 1, . . . ,d 
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and a positive real number p so that the following Agler decomposition holds: 

d 

p'-il^B)- S{5) {I ® B) Sir)* = J2 ^)[^ - ^'^^^fe] 

k=l 

for all B e C{K) and 5 = {8i, . . . , 5d), r = (xi, . . . , Td) in D„c,<i- 



One of the main results of ^28J is that the noncomniutative Schur- Agler class has 
a contractive Givone-Roesser realization. 

Theorem 6.6. (See ^^.) A given function F: D^cd -> «> /C, ^ ® /C) is 

in the strict noncommutative Schur- Agler class 5^°^, ^(i^, if and only if there 
exists a strictly contractive colligation matrix 



M = 



'A 


B 










C 


D 




U 




y 



for some Hilbert state space X — X\ 



Xd so that the evaluation of F at 



S = {di,...,5d) e 
F{S) = 

where 



.^d is given by 



Z{S) 



- Z{S)iA< 
Si 



. I^)y^Z{S){B (E) I,c) 



(6.16) 



' Sd 



Hence a function F: D„c,(i — * ^{1^ (§) /C, 3^ /C) is in the strict noncommutative 
H°° -space H!^''^(C{U , y)) if and only if there is a bounded linear operator 



such that 



so that F 



IS given as m 



A B 
C D 

A 

(I6l6l) . 





®i=iXk 




®i=iXk 




u 




y 



B 



< 1 for some p > 



liU and y are finite-dimensional Hilbert spaces, we may view SA'^^ ^^(W, y) and 
H^^°^{C{lA,y)) as matrices over the respective scalar-valued classes SA'^^d 



When this is the case, it is natural to define rational versions of SA°r 



and to consist of those functions in SA"^^ d (respectively, H^^ d) which 

the realization (|6.16p can be taken with the state spaces Xi, . . . ,Xd also finite- 
dimensional; we denote the rational versions of SA'^^ d ^^'^ ^^c d T^^-^nc d ^^id 
7?.iJ^'^, respectively. We remark that as a consequence of Theorem 11.1 in [57], this 
rationality assumption on a given function F in i?,'^'^ can be expressed intrinsically 
in terms of the finiteness of rank for a finite collection of Hankel matrices formed 
from the power-series coefficients Fa of F, i.e., the operators F^ G £{U,y) such 
that 

F{6) = V 



E 



In general, the embedding of a noncommutative integral domain into a skew- 
field is difficult (see e.g. [7S1 [52] )■ For the case of TZH°°'°, the embedding issue 
becomes tractable if we restrict to denominator functions D{S) € H"^ ° S ^0,(1/1) for 
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which -D(O) is invertible. If D is given in terms of a strictly contractive reahzation 
D{5) = D + C(/ - Z(J)A)-iZ(^)B (where A = A /k; and similarly for B, C 
and D), then D{8)^^ can be calculated, at least for ||Z(5)|| small enough, via the 
familiar cross-realization formula for the inverse: 

D{S)-^ = - D-iC(/ - Z(J)A^)-iZ((5)BD-i 

where A^ = A"" ® Ik. with = A-BD~^C. We define Q{nH^^'°^){C{U,y))o to 
be the smallest linear space of functions from some neighborhood of in D„c,(i (with 
respect to the Cartesian product operator- norm topology on D„c,d C C{IC)'^) to 
C{U, y) which is invariant under multiplication on the left by elements oiTZH^^''^{C{y)) 
and by inverses of elements of TZH'^^'°^{C{y)) having invertible value at 0, and invari- 
ant under multiplication on the right by the corresponding set of functions with lA in 
place of 3^. Note that the final subscript in the notation Q{TZH^^''^){£{U ,y))o is 
suggestive of the requirement that functions of this class are required to be analytic 
in a neighborhood of the origin G ©nc,d- 

Let us denote by TZO^^ ^{C{U, y)) the space of functions defined as follows: we 
say that the function G defined on a neighborhood of the origin in Bncd with values 
in £(W, 3^) is in the space TZO^^^ di^^^ ^)) '^^ ^ ^ realization of the form 

G{6) = D + C(/ - Z(J)A)-iZ((5)B 

for a colligation matrix M := ^] of the form M = M ® /yc where 



'A 


B 










C 


D 




U 




y 



for some finite-dimensional state-spaces Xi, . . . , Xd- Unlike the assumptions in the 
case of a realization for a Schur-Agler-class function in Theorem 16.61 there is no 
assumption that M be contractive or that A be stable. It is easily seen that 
Q{TZH^^''^{£{U ,y)))o is a subset of T?.©^^ ^(£(W, y)); whether these two spaces 
are the same or not we leave as an open question. We also note that the class 
TZO^^ di^i^^y)) tias an intrinsic characterization: F is in TZO^^ j;(£(W,3^)) if and 
only if some rescaled version F(5) = F(rS) (where rd = {rSi, . . . ,rdd) if ^ = 
(^1, . . . , Sd)) is in the rational noncommutative i7°°-class TZH^^''^{C{U , y)) for some 
r > and hence has the intrinsic characterization in terms of a completely positive 
Agler decomposition and finite-rankness of a finite collection of Hankel matrices as 
described above for the class TZH^ °{C{L(,y)). 

We may then pose the following control problems: 

Noncommutative polydisk internal-stabilization/iJ°°-control problem: We 

suppose that we are given finite-dimensional spaces W, U, Z, y and a block-matrix 
G = [^11 %\l ] in nOl^ di'^iW (BU,Z® y)). We seek to find a controller K in 
7?.0^^ ^(£(3^, W)) which solves the (1) internal stabilization problem, i.e. so that 
the closed-loop system is internally stable in the sense that all matrix entries of the 
block matrix 6(G, K) given by (|2.2[) are in "^77^'^, and which possibly also solves 
the (2) H°° -problem, i.e., in addition to internal stability, the closed- loop system 
has performance in the sense that T^w = Gn + Gi2K[I — G22K)~^G2i is in the 
rational strict noncommutative Schur-Agler class 7^5^° (W, Z). 

Even though our algebra of scalar plants T^-O^^, ^ is noncommutative, the param- 
eterization result Theorem 13.51 still goes through in the following form; we leave it 
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to the reader to check that the same algebra as used for the commutative case leads 
to the following noncommutative analogue. 

Theorem 6.7. Assume that G e TZO°^ ^{C{W ®U,Z Q) y)) is given and that 
G has at least one stabilizing controller K<,. Define = (I — 6*22^^*)^^, V* = 
K,{I -G22K,)-^, = {I-K,G22)~^ andV, = {I - K.,G22y^K,. Then the set 
of all stabilizing controllers K for G is given by either of the two formulas 

K ={V^ + Q)[U^ + G22Qy^ subject to (C/* + G22Q)(0) is mvertible, 

K = {U., + QG22y\V, + Q) subject to {U^ + QG22)iO) is mvertible, 

where in addition Q has the form Q = LAL where L and L are given by (j3.8p and 
A is a free stable parameter in H^^°^{C{y ®U,IA ® y)). Moreover, if Q = LAL 

with A stable, then {U^ + G22(9)(0) is invertible if and only if ([/* + QG22)(0) is 
invertible, and both formulas give rise to the same controller K . 

Given a transfer matrix G22 G TZO^^^ ^{C{U ,y)), we say that G22 has a stable 
double coprime factorization if there exist transfer matrices D{6), N(S), X{S), 
Y{5), D{6), N{S), X(d), and Y{S) of compatible sizes with stable matrix entries 
(i.e., with matrix entries in Ti-H'^^ "^) subject also to 

£1(0), 5(0), X{Q), X{Q) all invertible 
so that the noncommutative version of condition (13. 9p holds: 
G22(^) - D{5)-^N{8) = N{8)b-\8), 

' D{8) -N{8) 
-¥{8) Xi8) 

Then we leave it to the reader to check that the same algebra as used for the 
commutative case leads to the following noncommutative version of Theorem 13.111 

Theorem 6.8. Assume that G G TZO^^ ^ is stabilizable and that G22 admits a 
double coprime factorization (|6.17p . Then the set of all stabilizing controllers is 
given by 

K{8) = {Y{8)+D{8)A{8)){X{8) + N{S)A{8))-^ 

= iX{8) + Ai8)N{8))-\Y{8) + Ai8)D{8)), 

where A is a free stable parameter from H^/^{C(U,y) such that X{Q) - iV(0)A(0) 

is invertible and X{Q) + A(0)-/V(0) is invertible. 

Just as in the commutative case, consideration of the iJ°°-control problem for a 
given transfer matrix G G T^O^^ ^(£(W ® Z/^, Z © 3^)) after the change of the design 
parameter from the controller K to the free-stable parameter A in either of the two 
parameterizations of Theorems 16.71 and 16.81 leads to the following noncommutative 
version of the Model-Matching problem; we view this problem as a noncommutative 
version of a Sarason interpolation problem. 

Noncommutative-polydisk Sarason interpolation problem: Given matrices 
Ti, T2, T3 of compatible sizes over Ti-H'i^^"^, find a matrix A (of appropriate size) over 
'^^ncd that the matrix S = Ti + T2AT3 is in the strict rational noncommutative 
Schur-Agler class USA^^^dO^, Z). 

54 



Xi8) Ni8) 
Y{8) D{8) 



•-ny 





(6.17) 



While there has been some work on left-tangential Nevanlinna-Pick-type inter- 
polation for the noncommutative Schur-Agler class (see [H]), there does not seem 
to have been any work on a Conimutant Lifting theorem for this setup or on how 
to convert a Sarason problem as above to an interpolation problem as formulated 
in [22] . We leave this area to future work. 

6.3. Equivalence of state-space noncommutative LFT-model and non- 
commutative frequency-domain formulation. In order to make the connec- 
tions between the results in the previous two subsections, we consider functions 
as in Subsection 16.21 but we normalize the infinite dimensional Hilbert space K, to 
be i"^ and work with d -I- 1 variables 8 — [5i, . . . , 5d+i) in instead of d. As 

pointed out in Subsection 16.11 we may without loss of generality assume that the 
last variable 5d+i is fixed to be the shift operator S on P . 

The following is an improved analogue of Lemma 14.131 for the noncommutative 
setting. 

Theorem 6.9. Suppose that the matrix function W G TZO^^^ ^_^_i{C{U,y)) has a 
finite- dimensional realization 

W{S) = D + C(/ - Z{S)A)-^Z{S)B, 

where 

A = A (g) , B = B (g) , C = C (g) /|2 , D = D (g) , 

which is both nc-Hautus- detectable and nc-Hautus-observable. Then W is stable in 
the noncommutative frequency- domain sense (i.e., all matrix entries of W are in 
^ncd+i) '^^'^ on/y if W is stable in the state-space sense, i.e., the matrix A is 
nc-Hautus-stable. 

Proof. If the matrix A is nc-Hautus-stable, it is trivial that then all matrix entries 
of W are in We therefore assume that all matrix entries of W are in 

-^nc'd-i-i- I* remains to show that, under the assumption that {C, A} is nc-Hautus 
detectable and that {A, B} is nc-Hautus stabilizable, it follows that A is nc-Hautus 
stable. 

The first step is to observe the identity 



S,{S) ■■= 



I - Z{5)A 
C 



{I - Z{S)A)-^Z{d)B = 



Z{S)B 
W{d)--D 



(6.18) 



Since W{6) — D is in Hl^^''^^-^^{C{h{,y)) by assumption and trivially Z((5)B is in 
Hncld+ii^i^^^))^ it follows that SiiS) is in ff^;°^,(/:(iY, <Y © y)). By the de- 
tectability assumption and Proposition 16.11 it follow that there exists an operator 
L = L ® with L : y X so that A + LC is nc-Hautus-stable. Thus 

Fi(<5) = (/-Z(<5)(A + LC))-i [ / -Z{5)L] 

is in H^^'°^^{C{X ® y,X)). Note that Fi{6)Si{S) = (/ - Z{6)A)-'^Z{S)B. The 
fact that both Fi and are transfer-functions over -ff^'^^i implies that 52 (^) = 
(7 - Z(5)A)-iZ(<5)B is in H^^°,^,{C{U, X)). 
We next use the identity 

[Z(^) S2{5)]:^[Z{5) {I-Z{5)A)-^Z{5)B\ 

^Z{S){I - AZiS))-^ [I - AZ{S) B] . (6.19) 



55 



Now the nc-Hautus-stabilizability assumption and the second part of Proposition 
Oimply in a similar way that 83(6) = Z(^)(/- Z(^)A)-i is in H^^-°^^{C{X , X)) . 
Note that 53 in turn has the trivial realization 

D' + C'(/-Z(^)A')"^Z(^)B' 



On 



where [^;b;] ^ [a;^;] ^j^^ aj^d [^:^,] = [jl]. Thus {A' , B' ,C' , D') 
{A, I, I, 0) is trivially GR-controUable and GR-observable in the sense of [ 
the other hand, by Theorem 16.61 there exists a strictly contractive matrix [^// ^ ] 
so that 

SsiS) = r"C"{I - Z(^)A")"^Z(^)B" 
for some r < 00. Moreover, by the Kalman decomposition for noncommuta- 
tive GR-systems given in J27, we may assume without loss of generality that 
{A" , B" ,C" ,0) is GR-controUable and GR-observable. Then, by the main result 
of Alpay-Kaliuzhnyi-Verbovetskyi in [Hj, it is known that the function S{S) = 
'^aeJ^d ® ^" uniquely determines the formal power series S{z) — X^aej^^ SaZ°' . 
It now follows from the State-Space Similarity Theorem for noncommutative GR- 
systems in [27] that there is an invertible block diagonal similarity transform Q G 
C{X',X") so that 



'A 


I 




'A' B'' 




■q-i 0' 




'A" B"' 




Q 


0" 


I 







C 




/ 




C" 







/ 



In particular, A~Q ^A"Q where A" is a strict contraction and Q is a structured 
similarity from which it follows that A is also nc-Hautus-stable as wanted. □ 

We can now obtain the equivalence of the frequency-domain and state-space 
formulations of the internal stabilization problems for the case where the state- 
space internal stabilization problem is solvable. 

Theorem 6.10. Suppose that we are given a realization 



GiS) 



(/-Z(J)A)-iZ(^) [Bi B2] 



G2l(^) G22WJ [D21 J ^ [C2 

for an element G G TZO''^-^^ ^j^^{L{W (QlA, Z ®y)) such that the state-space internal 
stabilization problem has a solution. Suppose also that we are given a controller 
K G TZO'^^ j^j^^{C{y,b()) with state-space realization 

K{8) = T>K + Ck{I - ZK{S)AKr'ZK{S)BK. 

which is both nc-Hautus-stabilizable and nc-Hautus-detectable. Then the controller 
K {Aif, B/f , , Da'} solves the state-space internal stabilization problem as- 
sociated with {A, [Bi B2] ,[02]' [oil 0^]} '^^'^ ^^^y */-^(^) solves the non- 
commutative frequency- domain internal stabilization problem associated with 

r<(Ji\ — Gi2(<5) 

^^'^) ^ [G2i(5) G22(5) 

Proof. By Theorem l6.31 the assumption that that the state-space internal stabiliza- 
tion problem is solvable means that {C2, A} is nc-Hautus-detectable and {A, B2} 
is nc-Hautus-stabilizable. We shall use this form of the standing assumption. More- 
over, in this case, a given controller K ~ {Ak,Bk, Ck, Dif} solves the state-space 
internal stabilization problem if and only if K stabilizes G22 . 

Suppose now that K ~ {Ak,Bk,Ck,^k} solves the state-space internal 
stabilization problem, i.e., the state operator Ad in (|6.14p is nc-Hautus-stable. 
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Note that the 3x3 noncommutative transfer matrix &{G, K) has reahzation 
@{G,K) = De + Ce(/ - Ze(<5)Ae)-iZe ^Bq with Ze(^) = Z,i{8) as in (EUD 
where 



" Ae 


Be 




■ Ae 


Be ' 


. Ce 


De 




_ Ce 


De _ 



) 1^2 



with 

Ae = 

Ce = 



A + B2DKC2 

BkC2 

'Ci + D12DKC2 

DkC2 

C2 



BwCk 
Ak 

D12CK 
Ck 




B2 





D(- 



Bi + B2DKD2 

BkD21 

Di + D12DKD21 

DkD21 

D21 



B2DK 




Bk . 




D12 


D12D 


I 


Kk 





I 



K 



(6.20) 



Now observe that Ae is equal to Ac;, so that all nine transfer matrices in 0(G, K) 
have a realization with state operator Ae — A.ci nc-Hautus-stable. Hence all matrix 
entries of S{G,K) are in 

Suppose that K(S) with realization K ^ {Ak,^KiCk,^k} internally sta- 
bilizes G in the frequency-domain sense. This means that all nine transfer ma- 
trices in Q{G,K) are stable. In particular, the 2x2 transfer matrix W :— 
@{G22,K) - &{G22,K){0) is stable. From (jg:^ we read off that W has real- 
ization 

^(^)= [ C2 J (^-ZeWAe) . 

By Theorem 16. 9[ to show that Ac; = Ae is nc-Hautus-stable, it suffices to show 
that { [ °g^=^ ] , Ac; } is nc-Hautus-detectable and that { Ac; , [ ^^J" ] } is nc- 
Hautus-stabilizable. By using our assumption that {A/^ , B^f , Cif , D^f } is both 
nc-Hautus-detectable and nc-Hautus-stabilizable, one can now follow the argument 
in the proof of Theorem l4.9l to deduce that { [ ] , Ac/} is noncommutative 

detectable and that {Acj, 
needed. 



BaDj, 



} is noncommutative Hautus-stabilizable as 



□ 



We do not know as of this writing whether any given controller K in the space 
^^nc 1^^^ ^ nc-Hautus-detectable/stabilizable reahzation (see the dis- 

cussion in the Notes below). However, for the Model-Matching problem, inter- 
nal stabilizability in the frequency-domain sense means that all transfer matrices 
Ti,T2,T3 are stable (i.e., have all matrix entries in -ffj^'^^^) hence the stan- 
dard plant matrix G = [ ] has a stable realization. A given controller K 
solves the internal stabilization problem exactly when it is stable; thus we may 
work with realizations K ~ {Ak,^k,Ck,^k} with Ak nc-Hautus-stable, and 
hence a fortiori with both {Ck,Ak} nc-Hautus-detectable and {Ak,Bk} nc- 
Hautus-stabilizable. In this scenario Theorem 16.101 tells us that a controller K{S) 
solves the frequency-domain internal stabilization problem exactly when any stable 
realization K ~ {A/f , B/f , Cj<-, D/f } solves the state-space internal stabilization 
problem. Moreover, the frequency-domain performance measure matches with the 
state-space performance measure, namely: that the closed-loop transfer matrix 
Tzw = Gil + Gi2[I — i4rG'22)~"^^G2i be in the strict noncommutative Schur-Agler 
class iSyi°c ■2)- In this way we arrive at a solution of the noncommutative 
Sarason interpolation problem posed in Section [ 
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Theorem 6.11. Suppose that we are given a transfer matrix of the form G 
[t'Jo]^ Hnc°d+ii'^0^ ®U.Z® y)) with a realization 

(/-Z(,5)A)-iZ(J) [Bi B2] 



T,{8) 


T2iS) 




Dn D12 




Ci 









D21 


+ 


C2 



(so C2(/ - Z(<5)A)-iZ(J)B = for all 5) where 



A 

Ci 
C2 



Bi 
Dn 

D21 



B2 
D12 





A 
C2 



Bi 
D21 



B2 
D12 




00,0 

nc,d-\ 



as usual. Then there exists a K £ 
noncommutative Schur-Agler class ^5^1°^, ^ 
with 2? as in (|4.1ip . satisfying LMIs: 



so that Ti + T2KT^ is in the strict 
if and only if there exist X,Y £ V, 





No 




'AY A* - Y 
C\YA* 
BI 

'A*XA - X 
BIX A 



AYCl 
CiYCl - / 
^11 

A*XBi 
BlXBi - I 



Bi 
Dn 
-I 

D*n 
-I 









<o, r>o, 



<o, x>o, 



and the coupling condition 



X I 
I Y 



> 0. 



Here Nc and No are matrices chosen so that 

Nc is injective and Im Nc = Ker [Sj ^12] '^'^'^ 
No is injective and ImTVo — Ker [C2 i?2i] • 

6.4. Notes. 1. The equality of ^ a (A) with /2a (A) where A is as in (|6.8p appears in 
Paganini's thesis [108] : as mentioned in the Introduction, results of the same flavor 
have been given in [371 EH [Ml IM [US] ■ Ball-Groenewald-Malakorn ^ show how 
this result is closely related to the realization theory for the noncommutative Schur- 
Agler class obtained in [55]. There it is shown that /LtA(A) < JI^iA) — /2a (A), 
where ^^^(A) is a uniform version of /xa(A). The fact that /xa(A) = 71^ (A) is 
the content of Theorem B.3 in [lOSj . Paganini's analysis is carried out in the more 
general form required to obtain the result of Proposition 16. II 

The thesis of Paganini also includes some alternate versions of Proposition 16.11 
Specifically, rather than letting each Sj be an arbitrary operator on one may 
restrict to such operators which are causal (i.e., lower-triangular) and/or slowly 
time-varying in a precise quantitative sense. With any combination of these refined 
uncertainty structures in force, all the results developed in Section [S] continue to 
hold. With one or more of these modifications in force, it is more plausible to argue 
that the assumption made in Section 16.11 that the controller K has on-line access 
to the uncertainties di is physically realistic. 

The replacement of the condition /i(A) < 1 by /i(A) < 1 can be considered as a 
relaxation of the problem: while one really wants fJ.{A.) < 1, one is content to ana- 
lyze /2(A) < 1 since /2(A) is easier to compute. Necessary and sufiicient conditions 
for /2(A) < 1 then provide sufficient conditions for /x(A) < 1 (due to the general 
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inequality /i(A) < /i(A)). In the setting of the enhanced uncertainty structure dis- 
cussed in this section, by the discussion immediately preceding Proposition 16.11 we 
see in this case that the relaxation is exact in the sense that /i(A) < 1 is necessary 
as well as sufficient for /x(A) < 1. In Remark 1.2 of the paper of Megretsky-Treil 
|99j , it is shown how the /i-singular- value approach can be put in the following gen- 
eral framework involving quadratic constraints (called the S-procedure for obscure 
reasons). One is given quadratic functionals (Tqi fij ■ ■ • > cr^ defined on some set L 
and one wants to know when it is the case that 

aj{x) > for j = 1, . . . ,^ =^ CTo(x) < for a; £ i. (6.21) 

A computable sufficient condition (the relaxation) is the existence of nonnegative 
real numbers ti, . . . , (t, > for j = I, ...,£) so that 

e 

c^oix) + Tjaj{x) < for all x e L. (6.22) 

The main result of [HS] is that there is a particular case of this setting (where L 
is a linear shift-invariant subspace of vector- valued L^(0, cxd) (or more generally 
L^Qj,(0, cx))) and the quadratic constraints are shift-invariant) where the relaxation 
is again exact (i.e., where (|6.2ip and (|6.22p are equivalent); this result is closely 
related to Proposition 16.11 and the work of [108] . A nice survey of the S-procedure 
and its applications to a variety of other problems is the paper of Polik-Terlaky 

m- 

2. It is of interest to note that the type of noncommutative system theory devel- 
oped in this section (in particular, nc-detectability/stabilizability and nc-coprime 
representation as in (|6.17p ) has been used in the work of Beck 36] and Li-Paganini 
[89j in connection with model reduction for linear systems with LFT-modcUed struc- 
tured uncertainty. 

3. We note that Theorem 16.81 gives a Youla-Kucera-type parametrization for 
the set of stabilizing controllers for a given plant G £ T^O^^ ^(£(>V © -E © 3^)) 
under the assumption that G22 has a double coprime factorization. In connection 
with this result, we formulate a noncommutative analogue of the conjecture of Lin: 
If G e TZO°^ JC{W ®U,Z ® y)) is stabilizable, does it follow that G22 has a 
double- coprime factorization? If G22 has a realization 

G22iS) = C2{I - Z{5)A)-^Z[5)B2 

with ^] — [c f ] ® nc-Hautus stabilizable and nc-Hautus detectable, then 
one can adapt the state-space formulas for the classical case (see |1041[85] ) to arrive 
at state-space realization formulas for a double-coprime factorization of G22 ■ If it 
is the case that one can always find a nc-Hautus stabilizable/detectable realization 
for 022, it follows that G22 in fact always has a double-coprime factorization and 
hence the noncommutative Lin conjecture is answered in the affirmative. However, 
we do not know at this time whether nc-Hautus stabilizable/detectable realizations 
always exist for a given G22 G 'R-O^^ di^^^ ^))- From the results of [27], it is known 
that minimal i.e., controllable and observable realizations exist for a given G22. 
However, here controllable is in the sense that a certain finite collection of control 
operators be surjective and observable is in the sense that a certain finite collection 
of observation operators be injective. It is not known if this type of controllability is 
equivalent to nc-Hautus controllability , i.e., to the operator pencil [/ — Z{5)A B] 
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being surjective for all d G C{£'^)'^^^ (not just 6 in the nonconimutative polydisk 
Dnc.d)- Thus it is unknown if controllable implies nc-Hautus stabilizable in this 
context. Dually, we do not know if observable implies nc-Hautus detectable. 

4. Theorem 16.91 can be viewed as saying that, under a stabilizability/detect- 
ability hypothesis, any stable singularity of the noncommutativc function W must 
show up internally as a singularity in the resolvent (/ — Z((5)A)^^ of the state 
matrix A. A variant on this theme is the well known fact for the classical case 
that, under a controllability/observability assumption, any singularity (stable or 
not) of the rational matrix function W{X) = D + XC{I — XA)~^B necessarily must 
show up internally as a singularity in the resolvent (I — XA)~^ of the state matrix A. 
A version of this result for the noncommutativc case has now appeared in the paper 
of Kaliuzhnyi-Verbovetskyi-Vinnikov [82 ; however the notion of controllable and 
observable there is not quite the same as the notion of controllable and observable 
for non-commutative Givone-Roesser systems as given in |27j . 

5. Given a function S{z) = X^neZ'' ^nz" (where z = (zi, . . . , Zd) is the variable 

in the commutative polydisk D'' and we use the standard multivariable notation 
z" = z^'i . . . z^" if n = (m, . . . , nd) G Z^), we know from the resuhs of [SIEIES, that 
S has a contractive realization S{z) — D + C{I — Z{z)A)Z{z)B. In light of the work 
of [28], we see that any such contractive system matrix : (©fc^i-^fc © — > 

(®fc=i'^fc ® y) can also be used to define an element S of the noncommutativc 
Schur-Agler class SAnc,d{l^,y)- 

S{d) = D + C(/ - Z{6)A)-^Z{S)B 

where [col = [ c B ] ® ■ Thus a choice of contractive realization {A, B, C, D} 
for the commutative Schur-Agler-class function S can be viewed as a choice of 
noncommutativc lifting to a noncommutativc Schur-Agler-class function S(<5); the 
lifting property is that 

S{zl) = S{z) (g) If2 where zl = {zilp^ . . . , Zdip) G Dnc,d if z = (zi, . . . , z^) G B*^. 

While the realization for the commutative function is highly non-unique, the real- 
ization for the noncommutativc function is unique up to state-space similarity if 
arranged to be minimal (i.e., controllable and observable as in [27]). Philosophi- 
cally one can say that evaluation of the function on the commutative polydisk D'' 
does not give enough frequencies to detect the realization; enlarging the frequency 
domain (or points of evaluation) to the noncommutativc polydisk DJ^^ ^ does give 
enough frequencies to detect the realization in an essentially unique way. 
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